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TECHNICAL NOTE 2^09 

SUMMARY OF METHODS FOR CALCULATING DYNAMIC LATERAL 
STABILITY AND RESPONSE AND FOR ESTIMATING 
LATERAL STABILITY DERIVATIVES 
By John P. Campbell and Marion 0. McKinney 

SUMMARY 


A summary of methods for making dynamic lateral stability and 
response calculations and for estimating the aerodynamic stability 
derivatives required for use in these calculations is presented. The 
processes of performing calculations of the time histories of lateral 
motions, of the period and damping of these motions, and of the lateral 
stability boundaries are presented as a series of simple straightforward 
steps. Existing methods for estimating the stability derivatives are 
summarized and, in some cases, simple new empirical formulas are pre- 
sented. Reference is also made to reports presenting experimental data 
that should be useful in making estimates of the derivatives. Detailed 
estimation methods are presented for low-subsonic- speed conditions but 
only a brief discussion and a list of references are given for transonic- 
and supersonic-speed conditions. 


INTRODUCTION 


Dynamic lateral stability has not received widespread attention in 
the past because it has not generally been a serious problem in the 
design of airplanes. Consideration of dynamic lateral stability has 
recently become more important, however, because current design trends 
toward the use of low aspect ratio, sweepback, and higher wing loading 
have, in many cases, led to unsatisfactory dynamic lateral .stability. 
Airplane designers are therefore finding it necessary to make such calcu- 
lations in connection with the design and modification of airplanes. In 
many cases these calculations are difficult to perform for designers who 
have had no previous experience in theoretical stability work because 
most of the published theoretical analyses are not presented in a form 
that is especially suited to the computation of dynamic stability. The 
estimation of the stability derivatives required in dynamic stability 
calculations has also been found to be difficult in many cases. Although 
theoretical and experimental data on these derivatives have appeared in 
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numerous publications, no single publication has presented methods for 
estimating the derivatives for all types of airplanes. 

One approach to a presentation of methods of calculating stability 
and estimating stability derivatives in a form suitable for use by 
designers was made by Zimmerman in reference 1. Although this report 
has proved to be of valuable assistance to designers in making dynamic 
stability calculations, recent trends in airplane design have caused its 
usefulness to be seriously limited. For example, the equations of refer- 
ence 1 do not include the product-of-inertia terms which have been shown 
by recent studies to be very important in some cases. (See references 2 
and 3. ) Moreover, the calculation of the time histories of lateral 
motions, one type of calculation that has been the subject of increasing 
interest in the last few years (references 4 to 7), is not covered in 
reference 1. The methods of estimating stability derivatives presented 
in reference 1 are also limited because they apply only to airplanes 
having unswept wings with an aspect ratio of 6 operating at speeds at 
which compressibility effects are negligible. The purpose of the present 
paper is to extend the methods of reference 1 to include the methods of 
computation which are of current interest to designers and to include 
methods of estimating derivatives for configurations and flight conditions 
which are now being considered. 

This' paper summarizes and reduces to simple straightforward steps 
methods for computing the time histories of lateral motions, the period 
and damping of these motions, and the lateral stability boundaries. 
Existing methods of estimating stability derivatives for a variety of 
airplane configurations are summarized and, in some cases, simple new 
empirical formulas are presented. Reference is also made to reports 
presenting experimental data that should be useful in making estimates 
of these derivatives. 

SYMBOLS 

All forces and moments are referred to the stability system of axes 
which is defined in figure 1. The following definitions apply to the 
symbols except where they are otherwise defined: 

m 

S 

c 

b 

Yd 


mass of airplane, slugs 
wing area, square feet 
wing mean chord, feet (b/A) 
wing span, feet 

span of that part of wing that has tip dihedral, feet 
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l 


h 

w 


s 

X 

d 


Z H 


z 


A 

A 

X. 


r 

r T 

t 

v 




tail length (distance from center of pressure of vertical 
tail to center of gravity, measured parallel to longi- 
tudinal stability axis; values of Z must be calculated 
for each angle of attack), feet 

average fuselage height at wing root; feet 

average fuselage width at wing root, feet 

vertical distance of quarter chord of wing root chord from 
fuselage center line, positive downward, feet 

nondimensional time parameter based on span (Vt/b) 

longitudinal distance rearward from airplane center of 
gravity to wing aerodynamic center, feet 

longitudinal distance from leading edge of vertical tail 
chord to horizontal tail aerodynamic center, feet 
(see fig. 6) 

vertical distance from horizontal tail to base of vertical 
tail., feet (see fig. 6) 

height of center of pressure of vertical tail above longi- 
tudinal stability axis; values of z must be calculated 
for each angle of attack, feet 

aspect ratio 

sweepback of wing quarter-chord line, degrees 

taper ratio (Tip chord/Root chord); also, differential 
operator in Laplace transform 

dihedral angle, degrees (see sketch of fig. 9) 

dihedral angle of wing tip, degrees 

time, seconds 

airspeed, feet per second 

radius of gyration about principal longitudinal axis of 
inertia, feet 

radius of gyration about principal normal axis of inertia, 
feet 
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k X 

k Z 


K-v 


K9 


K 


•X 


K17 


s xz 


%Z ' b 2 


xz 


Ki = 


K xz 


K- 


X 


K- 


K 2 = 


XZ 


K„ 


7 

a 


P 

0 


radius of gyration about X axis, feet 
|| k X o 2 cos 2 Ti + k Zo 2 sin 2 T] j 

radius of gyration about Z axis, feet 




%l h 


7 l 

V’ 


v b 


kv /b 


A 


k z /b 


product- of- inertia factor lfk z ^ j sin cos rj j 


angle of attack of principal longitudinal axis of inertia, 
degrees (see fig. 2) 

angle of climb, degrees (see fig. 2) 

angle of attack of longitudinal body axis, degrees 
(see fig. 2) 

angle between principal longitudinal axis of inertia and 
longitudinal body axis, degrees (see fig. 2) 

air density, slugs per cubic foot 

angle of bank, radians 


\ 


J 
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♦ 

3 

P 

r 

to 

*0 

00 

(^)o 

(DlOo 

R 

I 

A,B,C,D,E 
p l> p 2> ••• p 7 


^*1 * ^2 * *3 > 

D 

P 

T l/2 


angle of yaw, radians 
angle of sideslip, radians 

rolling velocity, radians per second (d0/dt) 
yawing velocity, radians per second (dt/dt) 
initial angle of bank, radians 
initial angle of yaw, radians 
initial angle of sideslip, radians 
nondimensional initial rolling velocity (d0/da) 
nondimensional initial yawing velocity (d^/da) 

Routh’s discriminant or real part of complex root R + Ii 
imaginary part of complex root R + Ii 

coefficients of the characteristic biquadratic equation 
factors of the B, C, and D coefficients 

roots of characteristic biquadratic equation 

differential operator (d/dcr) 

period of the lateral oscillation, seconds 

time to damp to one -half amplitude, seconds 


T 

o 

N c 

Y c 


time conversion factor (m/pSV) 
nondimensional time factor ( t/r ) 
relative density factor (m/pSb) 
impressed rolling moment, foot-pounds 
impressed yawing moment, foot-pounds 
impressed lateral force, pounds 
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C Z f 


c nc 

C *c 

cl 

c D 
c i 
c n 

C Y 

q. 


impressed rolling -moment coefficient 
impressed yawing -moment coefficient 
impressed lateral-force coefficient 
lift coefficient (Lift/qS) 
drag coefficient (Drag/qS) 

rolling-moment coefficient (Rolling moment/qSb) 

yawing -moment coefficient (Yawing moment/qSb) 

lateral-force coefficient (Lateral force/qS) 

/I p 

dynamic pressure , pounds per square foot I 


Sc L 

Cl o * 5T 


K)a * 


5a 


Cd 0 - C » - S 

Sc 7 

Cz p = 5T 


Cd 

2 


Cn 0 


SC 

W 


n. 


SCy 

% = sr 


SC 2 

Cl P ~ ^ 
2V 


SC 


n 


'”p " 

°2V 


nA 


5Cy 

Cy P ” ^ 
. °2V 
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_ _ 5£i 

l-r -,rb 


& 


2V 


2n r ^rb 
°2V 

ac Y 

C *r =C7b 

°2V 


dCZr 

^Pr - dr 


_ 

ZP = 2K X 2 


^ C n Q 


n fl = 


P 2K Z 2 


C Yp 

^ = “2" 

Z = ^- 

P 4K X 2 




C JA 

_ J2_ 


4K Z 2 


y^ = 


% 

5T 


z = 


4K-V 


nr 


'n Y 


4Kn 


7 
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*< 

4 

II 


dCz 

l c = 2 

2% 


H c n c 

n c “ 2 

2K Z 2 


o 

II 

1 ° 
ro 

o 


i — i 

increment 

forces 

CVI 

increment 

lift 


in Cnp produced^ by- 
in Cn p produced by- 


lift and induced -drag 
drag not associated with 


H 

o 

Subscripts: 

wing 

fus 

tail 

design 

data 

exp 

V-tail 

e 


horizontal tail 
section lift curve slope 

wing 

fuselage 

used to designate vertical tail 

used to designate design under consideration 

used to designate design for which force-test data are 
available 

experimental 

V-tail 

effective 


H 


horizontal tail 
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CALCULATION OF LATERAL STABILITY 

AND RESPONSE 


Various types of calculations may be performed to indicate in some 
way the stability of an airplane or the response to gust disturbances 
and control manipulations* The calculations most commonly made are cal- 
culations of time histories of disturbed motions , period and damping of 
the free motions, and spiral and oscillatory stability boundaries (lines 
of neutral damping of the spiral mode and of the lateral oscillations). 
Step-by-step procedures for performing these types of calculations are 
explained in the text and derivations and additional pertinent material 
are presented in appendixes A to D. 

The period and damping calculations are the easiest of the three 
types to perform. For this reason, and because the dynamic lateral 
stability of airplanes is at present specified in the flying-qualities 
requirements in terms of the period and damping of the lateral oscilla- 
tion, period and damping calculations are probably the most commonly 
performed. 

Recent dynamic stability work has indicated, however, that the 
period and damping characteristics of the free motions of an airplane 
are not always a sufficient indication of whether the dynamic behavior 
of an airplane following various types of disturbances will be con- 
sidered satisfactory. For this reason the calculation of time histories 
of the motions of airplanes is becoming more common despite the fact 
that these calculations are fairly laborious. The increasing use of 
automatic computing machines has also made the calculation of motions 
more popular. 

For many years, calculations of stability boundaries were the type 
of calculation most commonly performed. In recent years, however, sta- 
bility boundaries have not been considered to give an adequate indica- 
tion of stability. Since boundaries are useful in some cases, however, 
(for example, for quick approximation of the effects of changes in 
dihedral and tail area) the methods of calculating the spiral and oscil- 
latory stability boundaries are described herein. Lines of constant 
period and damping of the lateral oscillation are related to stability 
boundaries (lines of neutral stability). In some cases these lines of 
constant period and damping may prove more useful than boundaries. Since 
no extensive use has been made of lines of constant period and damping, 
however, the methods of calculating these lines (presented in refer- 
ences 8 and 9) are not given in the present paper. 

The equations and methods of calculation presented in the present 
paper deal specifically with the inherent motions of airplanes for the 
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case of three degrees of freedom (roll, yaw, and sideslip) and linear 
stability derivatives. In order to perform similar calculations for 
cases involving additional degrees of freedom, nonlinear derivatives, 
or autopilots with time lag, special equations are required. The 
methods and equations for treating these cases are presented in refer- 
ences 10 to 18. Additional degrees of freedom for the case of free 
controls are treated in references 1 6 to 18 and for the case of fuel 
sloshing are treated in reference 10. The use of nonlinear derivatives 
in stability calculations is covered in reference 11. Methods of 
treating the effect of autopil9ts, including the effect of time lag in 
the autopilot are presented in references 12 to 15 and 19 . 

For some cases the effects of aerodynamic time lag are important. 
There are two different sources of such lag: (1) the time required for 

an aerodynamic impulse to travel from one component of the airplane to 
another (for example, the time required for a change in sidewash at the 
wing to reach the tail - a phenomenon commonly referred to as lag of 
sidewash); and (2) the time required for the growth and decay of the 
aerodynamic loads on the airplane components. For both of these cases 
the time-lag effects usually become increasingly important as the period 
of the lateral oscillation decreases. The effects of the first type of 
time lag can be accounted for in some cases by modification of the sta- 
bility derivatives. For example, the effect of the lag of sidewash on 
the derivative C nr is discussed subsequently under the section on 

"Estimation of Lateral Stability Derivatives". In many cases, how- 
ever, both types of time lag will require special stability equations. 

No general treatment of these cases has been published but an indi- 
cation of the method of treatment may be obtained from the treatments 
of autopilot lag in references 13 and 15 . 


CALCULATION OF PERIOD AND DAMPING 


As pointed out in references 1 and 2, the period and damping of the 
various modes of the lateral motion may be calculated from the roots of 
the characteristic equation 

AX 1 * - + BA. 3 + cx. 2 + DA. + E = 0 


by the equations 


P = 


2n 

I 


T 


and 



iogg 2 
R 


T 


0-693 


T 




R 
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where R represents a real root X or the real part of a complex root 
X = R ± Ii and I represents the imaginary part of a complex root. 
Negative values of T-^/2 represent the time required to double amplitude 

for unstable modes of the motion. 

The values of the coefficients A, B, C, D, and E may be 
obtained by the method given in steps 1, 2 , and 3 of the section on 
"Calculation of Motions". If the period and time to damp are to be cal- 
culated for a number of related cases, however, the values of the coef- 
ficients A, B, C, D, and E may be more conveniently calculated by 
a tabular procedure such as that shown as table I for making boundary 
calculations. 

Methods of determining the roots of the biquadratic characteristic 
equation are presented in appendix C . 


CALCULATION OF MOTIONS 


Calculation of the lateral motions of an airplane involves the 
integration of three simultaneous differential equations (see 
appendix A) to obtain a general solution in terms of the mass and 
aerodynamic parameters of the airplane. The general equations, once 
obtained, can then be used to obtain numerically the motions of any 
airplane in terms of the variation with time of the angles of bank, yaw, 
and sideslip or some function of these angles such as rolling or yawing 
velocity. Various methods, such as those given in references 20 to 22, 
are of course available for integrating the differential equations. 

Since the problems met in airplane dynamics are fairly complex, however, 
many of these methods are not suitable because of the difficulties of 
computation that arise. The method given in reference 4 (based on the 
Heaviside operational calculus) is satisfactory for calculating the 
forced motions following application of external forces or moments but, 
without modification, this method cannot be used to calculate the motions 
resulting from initial displacements in bank, yaw, or sideslip or from 
initial values of rolling Or yawing angular velocity. A solution based 
on the Laplace transformation is more satisfactory than that based on the 
Heaviside operational calculus because it permits direct calculation of 
the free motions following any initial condition, in addition to calcu- 
lation of the forced motions following application of external forces and 
moments. The application of the Laplace transformation to the calculation 
of lateral motions is outlined in appendix B. The material presented in 
this appendix is similar to the work presented in references 5 and 6 
except that the mass and aerodynamic stability derivatives have been com- 
bined as shown in appendix A to reduce the number of arithmetical and 
algebraic processes required in numerical solutions. 
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The process of calculating the motions is presented as a series of 
simple though lengthy arithmetical and algebraic steps so that an under- 
standing of the calculus involved in solving the differential equations 
is not required. The method as shown is suitable for calculating the 
motions as variations of 0, 3, p, and r with time for the case 

of the free motions following initial angular displacements (0 O , ty Q , 
and 3 q ) and angular velocities (D0) o and (lty) 0 and for the case of 
the forced motions resulting from constant impressed forces and moments 
(L c , N c , and Y c ). These are the cases for which motions are usually 
calculated. It is also possible to calculate the motions resulting from 
impressed forces and moments which are arbitrary functions of time by 
the methods explained in references 6 and 7« 


Motions Resulting from Initial Angular Displacements and Angular 

Velocities and from Constant Impressed Forces and Moments 

The six steps involved in obtaining a specific solution for the 
lateral motions of an airplane are : 

Step 1: Determine values of the following parameters: 

(a) Mass characteristics: 

k Xo> k Z 0 > 11 , and p 

(t>) Geometric characteristics: 

S and b 

(c) Flight conditions: 

V, C L , and y 

(d) Aerodynamic stability derivatives: 


Cl y C n y Cy y Cl y C y Cy y 

L $ } n 3* x 3* L -p x p' 



C n , and 

iip 



The methods of determining the values of the aerodynamic stability 
derivatives are given in subsequent sections of this paper. 

In cases where impressed forces and moments are used as disturbances, 
determine the values of the factors 

c v V % 

that are appropriate to the particular problem. 
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Step 2: From the known factors, evaluate the following parameters 

which are the stability derivatives in the form in which they are used 
in the calculation of motions! 


K XZ 

*1 K 2 
K X 


K xz 

2 '^ 


m 

T “ pSV 

* * ^ 

C 7 

np ' 2K Z 2 

Cll P 

yp = 

= V 

c h 

II 

c n p 

y P = 

c-i 

4 

11 

ro 

% 

” r ' >*# 

S 

y r = 


Also, when impressed forces and moments axe used, evaluate 


Z c = c z 
2% 2 ° 


2 C nc 


y c 2 C Y c 


2 2 

The values of , and Kyr? can determined from the following 

expressions 

2 2 2 2 2 
Kx = Kx 0 cos t) + KZq sin t) 


= K* ^cos^T| + Ky ^sin^T) 
K XZ = ( K Z 0 2 - % 0 2 ) sin i\ cos *1 


where 


tvy = , 

A 0 b 


YL<7 = 

■ o b 

Step 3: Solve for the values of the appropriate ones of the 

following coefficients from equations (l) to ( 4 ): 


Ik 
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In all cases solve for the values of A, B, C , D, and E: 
A = 1 - K x K 2 

B = Pi - Ay p 


where 


C = 


D = 


E = 


P iyp + p 2 + p 5yp + p 6y r - p 6 

7 


cl cl 

p 5 T + p 6 T ttul 7 + F 


c l Cl 

P 3 ~2 + p 4 ~2 tan V 


J 


P 1 z p " ^ + + 


P 2 = l jPr - l r*p 


p 3 _ ^p 1 ^ “ z r n p 
p 4 = z p n P " ^“p 


p 5 “ Kl^p ' z p 


p 6 = K 2 z p - n p 

p 7 = - p 2yp + p 3y p + p 4y r - p 4 


(1) 


The quantities to P7 are factors of the coefficients B, C, T), 

and E which are combinations of terms that occur frequently in calcu- 
lations of motions resulting from initial angular displacements and 
velocities and which are consequently grouped together for convenience. 
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Calculate the values of a^ a 1# . . . a^ when solving for the 
angle of bank 0 or the rolling velocity p: 


00 = 0oA 


aj. = 0o B + ^)o A 


&2 ~ 0 O C “Po p 5 + ( D ^)o(" A yp + Kfi J r ~ n r) “ 

(D^) 0 (K 1 n r - l r ) + l c - n c K! 

a 3 = 0 O ( P 6 ~t tan 7 + p t) ~ ^o p 5 ^ tan r - 3 0 p 3 + (D0) o ( p 6y r " 
p 6 - *2 l ryp + % y p) + (W) o(' P 5 y r + P 5 + ^Vp ‘ Z ^p) ' 
^(“r + yp) + n c( K l y P + l r) - y c P 5 

p 

a 4 = (0 o p 4 - *o p 3 + ( p 0) o p 6 - (^) 0 p 5)t tan 7 + 

*c( n p - n P y r + "r^ + n c( z 3 y r " Z p " Z r y p) " y c P 3 


V (2) 


= (" Z c n P + Vp)^ tan 7 


J 
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Calculate the values of bg , b^ , . . . , b^ when solving for 
angle of yaw ^ or the yawing velocity r: 

b 0 = V "N ' 

b x = V 0 B + (Dt) 0 A . 

b 2 - V ' >o P 6 - Wok 1 ! ' ”p) + (rt >o(-% * K 1“P * ^ - 
Z c k 2 + n c 

b 3 = ^o p 6 T + ^o( p 5 T + p t) “ p o p 4 + ( p 0) o (- p 6y p + *2^ - 

Vp) + (Dt) o(v P - k iVp + Vp) + l c(Vp + ”p) * > (3) 

n c( l p + y p) " y c p 6 

b 4 = [;0o p 4 + V3 - (°^o p 6 + - 

°P y p) + n c (Vp ■ Vi) • y c p 4 

b 5 = ( Z c n p ~ n c z p)ir 
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Calculate the values* of c n , c 


0' u l> • 


angle of sideslip (3: 


, "when solving for the 

.... \ . ■ ' ■ \ i 4 *• 


c 0 - Pq A 


"\ 


C 1 = 0o A ~t + V ^ tan 7 + 3 0 P! + (D0) o Ay p - (D*) 0 A(y r - 1) + y c A 


C L 


c l 


c 2 - Vl T * Vi T tm r + p o p 2 * W>o£ir - *aVr + Vp + 

Vr " "p + ( K 2 ! r - n r) y J * <D *>o[ A IT 7 + K lVr ‘ 

K 1 ”P * Vr + l p - (*l”r - l r^p] + ! = (‘Vr + K 2 + *p) + 

n c^r ^ ^1^^ + ^c p l 


Cl Cl , / Cl Cl 

00^2 2" + ^o P 2 2~ tan 7 + ( D 0)oi _K 2 Z p 2 ~ tan 7 + ~ tan 7 + 

Cl Cj\ f ° Cl 

K 2 z r ~ - ^ ~) + (Dt)o(Kin p -g- tan 7 - l p tt tan 7 - 

Cf CjN / Cl Cl \ 

K l% 2" + Z r in + z cf n p yr - °p - ^p + 2" - K 2 ^ tan 7 J + 

c l \ 

+ tan 7 I + y c P 2 


>( 4 ) 


M-Vr + Z P + z r y P - % 


/ c L c L \ / c L . c L \ 

c 4 = Z C^P T tan 7 ' °r 27 + n c (Jr V - r p T tan y 

Step Solve for the roots A.g, A.^, and A,^ of the 

biquadratic equation 

AA.^ + BA,3 + CA. 2 + DA, + E = 0 


s . 


(5) 
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where the values of the coefficients A, B, . . etc. were given by 
the solution of equations (l). Methods of determining the roots of the 
biquadratic equation sire given in appendix C. 

Step 5: Use the coefficients obtained from equations (l) to (4) 
and the roots of equation (5) to solve for the following coefficients: 


Calculate the values of the factors A^, A^, . . . , Ag when 

solving for the angle of bank 0 or the rolling velocity p: 


SqXi^ + 


+ 

a^X ^ 2 + 

a 4^l + 


6 AX X 5 

4 

+ 5 BX X 

+ kCX-^ 

2 

+ 3DX-L 

+ 2EX X 


SqX 2 *^ H 

h ^ 


0 

- a 3 X 2 h 

■ ajjXg + 

a 5 

6 AXr/ 

5 + 5 BX 2 l 

y + 4CX2 : 

5 + 3 DX 2 £ 

' + 2EX 2 

. V-3 5 + 

+ 

a 2^3^ + 

a 3 l 3 2 + 

V 3 + 


6 ax 3 5 

4 

+ 5 bx 3 

+ 4cx 3 3 

+ 3 DX 3 2 

+ 2**3 


a 0^4^ + 

a lH^ + 

a.2^1^ 

a 3 X 4 2 + 

ajjXjj. + 


6 ax ^ 5 

4 

+ 5 BX 4 

+ 4 CX 4 3 

+ 3DX 4 2 

+ 2 EX 4 





( 6 ) 
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Calculate the values of the factors B 2 , . . . , Bg when 

solving for the angle of yaw ^ or the yawing velocity r: 


Bn = 


+ h^X^^ + B 2 X^3 + b^X^^ ^*4^1 + h^ 

6AX3 5 + 56X3^ + ^CX!^ + 3DX3 2 + 2EX x 


Bo = 


^0^2^ "^2^2^ t b^X^ ^4^*2 t 

6ax 2 5 + 5BX 2 ^ + 4cx 2 3 + 3 DI 2 2 + ^^2 


B 3 = 


R L. 3 p 

b Q X 3 + b-^X 3 ^ b 2 X 3 + ^2^3 **" 


b|^X 3 


+ be 


6AX 3 5 + 5BX 3 4 + kcx^ 5 + 3BX 3 ^ + 2EX 3 




( 7 ) 


B4 




6AXlp + 5B Xl^ + 


2 

b 3 X]^ + b^Xj^ + b^ 

+ 3DX4 + 2EXj| 
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Calculate -the values ,of the factors. Cj_ , , C 2; , • 
solving for the angle of sideslip. . 0 : . • 


when 


C 1 = 


c 0 Xi5 + c^X^^ - c 2 x 1 3 + c^X^^ + c^Xj- 

6AX-l 5 + ^BX-j^ + l4CX x 3 + 3DX-^ 2 '+ 2 EX-l 




C 2 = 


C 0^2^ c 1^2 C 2^2^ ci 3 X 2 c j^ X 2 

6AX2 5 + 5BX2 4 + 4CX 2 3 + 3DX 2 2 + 2EX 2 


Co = 


; 0 X 3 ^ ^ C1X3 1 * + 02X^3 


+ C 3 X 3 d + c 


4 X 3 


3 6AX35 + 5BX3 1 * + Ucx 3 3 + 3DX3 2 + 2EX- 


> 


( 8 ) 


CqX^ 5 + C X X^ + c^ 3 + C3X^ 2 + c^X^ 

• 6a X4 5 + 5BX^ 4 + 4 CX^ 3 + 3DX]^ 2 + 2EX^ 


C 5 ~T 


J 


If equation (5) has conjugate complex roots, the values of the 
coefficients (equations ( 6 ) to ( 8 )) corresponding to these roots will 
he conjugate complex. In order to facilitate treatment of this case it 
is convenient to establish some special notation. This special notation 
is explained in appendix D. - 
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Step 6: ' The equations' of motion are written in different' form 
depending upon the roots of equation ( 5 ) ■ If the characteristic equa- 
tion has four real roots Xp, X2, X3, and XI*, the general form of 

the equations of motion is used, as follows: 


0 

* 

P 

P 


0X 4 


UM u U Ao 

Ai e x + Age + A^e + A^e ^ + AcO + Ag 
aXi crXn aXn aXh 

B-^6 4 * Bg6 4 ~ B36 ^ B^6 4 * B zzp 4 - Bg 

' aXn cXo aXo aXh 

C-]^6 + Cg6 4* C^e 3 4. Cjj6 4 - 

\( aXi crXp ctXo crXh \ 

dA^X-j^e + AgXge 4- ^3^3^ ^ A]^X]^6 + Ac/ 

1/ aXi aXp ctXo crXj. \ 

r\BlXi e + BgXge + ^3X36 ^ + Bj^Xj^ 6 + Bty 






j 


(9) 


If, as is generally the case, equation ( 5 ) 
real roots (R + Ii, R - Ii, X3, and X4), 
expressed as 


has two complex roots and' two 
the equations of motion may he 


, 0R . . 0X-, 0X1, 

0 = K^e cos( 0I + <Ofi) + A^e + A4e 4 + A^a + A^ 


A 


= Kge 0 ^ cos (0I + d>g) + B^e 0 ^ + B4e ff ^4-+. B^o + B^ 

0 = I^e 0 ^ cos (ol + coq) + C^e 0 ^ + C4e a ^ + 

p = ^R 2 + ^ e °^ cos ^®^ + to^ + 4 an"^ f) + 

A 3 X 3 e<J ^3 + A^\^e a ^ + A^J 

1 C X 2 2 oR f -lrt 

r = f KgYR + I e coslol + (Dg + tan gy + 

0X, 0X4 1 

63X36 ^ 4 * B||X]|6 4 - BeJ 


> do) 
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where 


K A = 2fH A 2 + I A 2 


K B = 2^/ + I B 2 


^ = W + X C 2 


tan- 1 ^ 

% 


+ - 1 Ib 
tan =- 


, -1 ic 

tan =- 

“c 


and R a and I A are defined in appendix D. 

If there are four complex roots (R + Ii, R - Ii, R* + I'i, and 
R ' - I ' i ) , the equations are 

0 = K A e° R cos(aI + a> A ) + K A 'e aR cos (cxl' + a> A ') + A^cr + Ag 

ijr = K B e aR cos(aI + o^) + K B 'e aR cos (ol* + o>g') + B^a + Bg 

4 

3 = K(-.e aR cos (al + a^) +'K^'e ffR cos(crI' + 

p = iL^ 2 + I 2 e aR costal + a) A + tan -1 | J + + I 


K A'p |2 + 1,2 e ° R cosfaI ! + co A ' + tan 


- 1 iij 

R'/J 


r = yJ^B^R 2 + ^ e° R cos + ‘“B + tan 1 g') + B^ + 
K b '|/r' 2 + I' 2 e aR ' cos^crl' + oig' + tan -1 |rj| 
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where 



cd a 1 = tan 


V 

V 


= tan 


= tan 


■1 ±A 
H a 


•1 IB 


■1 ic 

Rp 


(Ha) 




The coefficients K A> k b> Kq, ca A , o>g, and cjoq are defined in equa- 
tions (10a) and R A , I A , R A *, and I A * are defined in appendix D. 


Solve the appropriate ones of these equations of motion (equa- 
tions (9)> (10), or (11)) by substituting values of the nondimensional 
time factor o in the equations and solving for 0 , \|r , (3, py or r. 


Motions Resulting from Arbitrary Disturbances 

The motions resulting from arbitrary forcing functions can be 
obtained from the motions resulting from constant impressed forces and 
moments by the methods explained in references 6 and 7* 

A very useful method of obtaining the motion resulting from various 
abrupt gust and control disturbances is given by Jones in reference 7. 

In this paper it is pointed out that, although the component motions of 
an airplane must be calculated simultaneously (that is, by simultaneous 
differential equations), the effects of component disturbances may by 
the principle of superposition be calculated separately and later added 
in any desired proportion. Thus, if a given rolling moment causes a 
20° bank in 1 second and if a given yawing moment causes a 5° bank in 
1 second, the combined effect of both acting simultaneously will be a 
25° bank in 1 second. Jones also points out a somewhat similar fact 
with regard to the effects of disturbances that are not applied simul- 
taneously. This fact is that, if a given disturbance which arises at 
the time t = 0 is later augmented, the effect of the increment of 
disturbance will run its course independently of the effect of the 
original disturbance. For example, in a problem involving the correc- 
tion for a gust disturbance by a manipulation of the control, the motion 
produced by the gust disturbance can be calculated independently and the 
motion caused by the assumed corrective control manipulation can be 
added to it at any desired point. This example is illustrated graphi- 
cally in figure 3. 
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The principle of superposition may he applied analytically as veil 
as graphically. The analytical application which makes use of Carson's 
integral or Duhamel's integral is described in references 7 and 23. 

This method is useful for calculating the motions resulting from 
impressed forces and moments which are arbitrary functions of time. By 
application of these methods, the solutions for constant impressed 
forces and moments can be used to obtain new solutions for any arbitrary 
variation of impressed forces and moments with time which can be 
expressed by a mathematical formula. Some simple variations of 
impressed forces and moments with time and their Laplace transforms are 
given in reference 6. The transforms for any other function for which 
transforms have been worked out may be found in tables of Laplace 
transforms . 


CALCULATION OF STABILITY BOUNDARIES 
Oscillatory Stability Boundaries 


As pointed out in the preceding section of this report, the degree 
of stability of the uncontrolled motions of an airpiane is indicated by 
roots of the characteristic equation 


AA,^ + BX 3 + c\ 2 + DX + E = 0 


For stability the real roots or the real part of the complex roots of 
the characteristic equation must be negative., A useful discriminant 
for determining some of the characteristics of the roots in stability 
work is Routh T s discriminant R (R = BCD - AD^ - B^E). The use of this 
discriminant in dynamic stability analyses has been pointed out in many 
reports, for example, references 1, 2, 3, 21, and 2k. Routh has 

shown (reference 20) that, if R and the coefficient E are finite, 
the necessary and sufficient conditions that the real roots and the 
real parts of the complex roots should be negative are that every coef- 
ficient of the biquadratic and also R should have the same sign. 

Routh also showed that when R = 0 and B and D have the same sign 
there are a pair of complex roots with the real parts zero. Since the 
value of the real part of a complex root indicates the stability of an 
oscillatory mode of the motion of an airplane, the lateral oscillation 
is' neutrally stable when R = 0 and the coefficients B and D have 
the same sign. Oscillatory stability boundaries can be determined, 
therefore, by solving the equation R =■ 0 and checking to determine 
whether the signs of B and D are the same . 


Since two of the most important stability derivatives affecting 
lateral stability are the directional stability derivative C n ^ and 
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the effective dihedral derivative Cjp, boundaries for neutral oscil- 
latory stability are usually calculated as a function of these two 
derivatives as illustrated in figure These calculations are gener- 
ally carried out by the method shown in table I. This table contains a 
numerical example and step-by-step instructions for using the table. 

The results of this numerical example are plotted in figure 4 . The 
procedure illustrated in table I is first to assume values of the inde- 
pendent variable * Cnp to cover the range for which the boundary is 
required. The values of all the other mass and aerodynamic stability 
derivatives except Cjp are then estimated. The value of Cnp is 
generally assumed to have been varied by varying the size of the verti- 
cal tail and consequently the tail contribution to each of the other 
stability derivatives varies, as Cnp .is varied. The values of the 
coefficients A, B, C, D, and E and then R are calculated as 
functions of Zp! 


Z 


P " 


JL 

4k x 2 



The values of Zp corresponding to the assumed values of C n g for the 

condition of neutral oscillatory stability are next obtained by solving 
the expression R = 0 which is a quadratic in Zp that is of the form 

p 

u 1^3 + v l z 3 + W]_ = 0 

Finally, the values of corresponding to the assumed values- of C n p 

are obtained from the values of Zp. 


The values of Zp which satisfy the expression R = 0 must be 
checked to determine whether they satisfy the other condition for 
neutral oscillatory stability - that the sign of the coefficients B 
and D must be the same. This check can be performed readily by sub- 
stituting the values of Zp which satisfy R = 0 into the expression 
for D which is a linear equation of the form 


D = u 2 Zp + v 2 


Thus, the sign of D 
any given value of C 

three predominant terms of 


is determined. The sign of B is a constant for 
np and is almost invariably positive since the 


B contain the derivatives C 


Z -n ; 


and Cyp which in all practical cases contribute a positive increment 
to the value of B. 
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Since tvo values of C z p satisfy the condition R = 0 for each 
value of Cnp> the R = 0 curve has two branches. As pointed out in 

reference 24, one of the branches of the R = 0 curve generally repre- 
sents an oscillatory stability boundary and the other branch represents 
a line of numerically equal real roots with opposite signs. (See 
fig. 4.) If neither of the values of which satisfy the expres- 

sion R = 0 for a particular value of C n p is found to represent a 

point of neutral oscillatory stability, the lateral motion has no oscil- 
latory mode for that value of C n p. If both of the values of Cjp 

which satisfy the expression R = 0 are found to represent points of 
neutral oscillatory stability, the lateral motion has two oscillatory 
modes. In this case, since the boundary D = 0 represents the line of 
infinite period, the branch of the R = 0 boundary which lies close to 
the D = 0 boundary is usually the boundary for neutral stability of 
the longer period of the two oscillatory modes. A detailed discussion 
of the significance of the stability boundaries and the regions formed 
by these boundaries is given in reference 24. 

In calculating stability boundaries for a specific airplane a com- 
plete solution such as that explained in the preceding paragraphs should 
be made. For general studies of stability, however, approximate oscil- 
latory stability boundaries may be calculated much more simply by the 
methods shown in reference 24. 

As pointed out previously, methods of calculating lines of constant 
period and damping of the lateral oscillation are presented in refer- 
ences 8 and 9* 


Spiral Stability Boundaries 


Spiral stability boundaries, like oscillatory stability boundaries, 
are usually determined as a function of the directional stability deriva 
tive C n p and the effective dihedral derivative as illustrated in 

figure 4. As pointed out in reference 1, neutral spiral stability 
occurs when the E coefficient of the characteristic equation is zero 
(E = 0). A spiral stability boundary can be easily obtained from this 
relation. If expressions for E (in terms of Zp) corresponding to 
several values of C n p have already been obtained in the process of 


calculating an oscillatory stability boundary, the equations formed by 
setting these expressions for E equal to zero can be solved for the 
Zp ^and hence C z p) corresponding to the assumed values 

If the values of E have not already been obtained* in the 


values of 
of C 


nr 


process of calculating an oscillatory stability boundary, a spiral 
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stability boundary for the level-flight condition (7=0) can be cal- 
culated simply from the equation 

c *3 - Cn P (12) 

Values of Cnp are assumed within the range for which the boundary is 
required. The values of C2 r and C nr corresponding to each value of 
C n p are then determined. The tail contributions to these derivatives 
generally vary with Cn^ since C n ^ is usually assumed to be varied 
by changing the size of the vertical tail. 

ESTIMATION OF LATERAL STABILITY 

DERIVATIVES 

GENERAL REMARKS 


Methods of estimating the lateral stability derivatives have been 
presented in numerous publications but no single report has contained 
information for estimating the contribution of all principal airplane 
components to all the derivatives for airplanes having any sweep angle 
or aspect ratio. In the present paper , an approach to such a presenta- 
tion is made by the coordination of and reference to existing estima- 
tion methods, by reference to publications containing data which should 
be useful in making estimates, and by the suggestion in some cases of 
simple new empirical formulas. Detailed estimation methods are pre- 
sented for low-subsonic-speed conditions but only a brief discussion and 
a list of references are given for transonic- and supersonic -speed con- 
ditions. In general, the estimation methods presented should be expected 
to yield only fairly accurate values suitable for making first approxi- 
mations of dynamic stability. This limitation applies especially to 
the cases in which the derivatives are based completely on theory. 

For convenience, the references that should be useful in estimating 
the stability derivatives are presented in table II. The references 
are grouped according to the speed range covered (subsonic or super- 
sonic) and according to the derivatives presented in each report. The 
references for the subsonic case (references 1 and 25 to 9*0 are further 
divided into two groups - one including reports which contain estimation 
methods and the other including reports which contain experimental data 
that should be useful in making estimates of derivatives. The 
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references for the supersonic case (references 95 to 115 ) are sub- 
divided according to wing plan form. 

The following sections covering the estimation of the nine sta- 
bility derivatives are divided into three groups according to the type 
of derivative - sideslip derivatives (Cy^, C^p, Cjpj, r °ll in 6 deriva- 
tives (Clip, Cjp, c Ypj, and yawing derivatives (c nr , C j r , Cy r ] . The 
derivatives Cyp and Cy r have usually been neglected in making 

dynamic lateral stability calculations because theory indicated that for 
unswept wings Cyp and Cy r were zero. Recent experimental, data, 

however, have indicated that both swept and unswept wings produce meas- 
urable values of these derivatives (references 25, 59, and 86). Since 
the vertical tail contributes to Cy p and Cy r , it appears desirable 

to estimate these derivatives and to use them in the calculations of 
stability unless it is established that for the case in question the 
effects of Cyp and Cy r on stability are negligible. For these two 

derivatives, only the effect of the wing and vertical tail need to be 
considered . 

The methods of estimating the rolling and yawing derivatives pre- 
sented herein were obtained from theoretical treatments based on the 
assumption of steady rolling and yawing and from experimental data 
obtained principally from tests made under conditions of steady rolling 
and yawing. The only information that applies directly to the oscil- 
latory case is a limited amount of data on C nr obtained by oscillation 

techniques. When calculations are made in which the oscillatory mode 
is the subject of interest, some consideration should be given to cor- 
recting the derivatives based on steady rolling or yawing to account 
for differences in the derivatives that are likely to exist as a result 
of differences between the oscillatory motion and the steady rolling 
and yawing motion. For example, the data of reference 82 have indi- 
cated that, for flap-extended or power-on conditions, fairly large dif- 
ferences might exist between the values of the tail contribution to Cn r 
for the steady yawing and yawing oscillation cases. At present little 
information is available for correcting the values of C nr for the 

steady yawing case to apply to the oscillatory case and, unfortunately, 
little or no information is available for correcting the other stability 
derivatives. 

Since most wind-tunnel force-test data that are likely to be used 
in making estimates of the stability derivatives are probably for much 
lower Reynolds numbers than those for the full-scale airplane, some 
adjustments to the data are usually required to account for the dif- 
ferences in Reynolds number. The effects of Reynolds number should be 
considered in the cases of all the derivatives, especially those which 
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are estimated by methods that involve the use of force -test data. 

Methods of correcting for Reynolds number effects for some of the 
•derivatives are discussed in the following sections which cover the 
estimation procedures. In the cases where the Reynolds number effects 
are not discussed, it can be assumed that any abrupt variation in the 
derivatives near the stall for low-scale data will also be present for 
the full-scale airplane but will probably occur at a higher lift coef- 
ficient because of the higher maximum lift coefficient of the airplane. 

An indication of the lift-coefficient range over which the theory may 
not be expected to give reliable values of stability derivatives for the 
full-scale airplane can be obtained from large-scale drag data. The 
analysis of reference 86 indicates that the variation of the derivatives 
with lift coefficient is different from the theoretical variation at 
lift coefficients above that at which the drag due to lift increases 
abruptly from the ideal value Cl^/itA. 

The effects of Mach number and power are not treated in the sections 
on the individual derivatives but are discussed briefly in separate 
sections. A detailed treatment of these effects, including design 
formulas and charts, was considered beyond the scope of this paper. 


THE SIDESLIP DERIVATIVES Cy p , C n ^, C 


No satisfactory purely theoretical methods have yet been developed 
for obtaining accurate estimates of the sideslip derivatives C n p, 

and for a complete airplane, primarily because of large inter- 

ference effects between the various airplane components and because of 
large, and often unpredictable, variations of the derivatives with angle 
of attack. Fortunately, these derivatives can be obtained from conven- 
tional wind-tunnel force -test data. Such experimental, data are essential 
to the accurate determination of sideslip derivatives. It is, of course, 
highly desirable to have force-test data for the exact airplane design 
under consideration, but reasonably accurate estimates can usually be 
made by correcting the force-test data for a generally similar design. 

The methods of correcting the force-test data on a similar design for 
use in the case under consideration are covered in the following sec- 
tions. In the formulas presented, the subscript word "design" is used 
to designate the design under consideration and the subscript word "data" 
is used to designate the similar design for which force-test data are 
available . 

Force -test data should be used to determine the effect on the side- 
slip derivatives of such airplane components as leading-edge high-lift 
devices, stall-control devices, trailing-edge flaps, nacelles, external 
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stores, canopies, and dorsal and ventral fins. The effect of leading- 
edge high-lift devices is usually merely to extend to a higher lift 
coefficient the same variation of the derivative with lift coefficient 
as for the plain wing. Trailing-edge flaps often have large effects on 
the contributions of both the wing and the vertical tail to the sideslip 
derivatives (references 39 and 69); and since these effects are not 
easily estimated, it appears that in these cases use of force-test data 
is essential. The addition of nacelles and external stores generally 
has been found to decrease the directional, stability factor C np 

slightly. The results of a limited amount of research to determine the 
effect on the sideslip derivatives of the size and shape of canopies 
has been reported in references 48 and 73 but these results are inade- 
quate for making accurate predictions of the effects of canopies. The 
effects on the sideslip derivatives of dorsal and ventral fins are 
usually small at the small and moderate angles of yaw that are generally 
considered in stability calculations. (See references 47 and 71 •) 



In estimates of the lateral force due to sideslip derivative Cyp, 

force-test data for the design under consideration should be used when- 
ever possible. If such data are not available, data for a similar 
design can be used and corrected as follows: 

Wing-fuselage . - Since the wing-fuselage contribution to Gyp is 

usually relatively small compared with that of the vertical tail, great 
accuracy is not required in estimating this factor. This contribution 
may be estimated as follows: 

(1) Wing: If the wings of the two designs are generally similar 

the difference in Cyn can be considered negligible and no correc- 

ting 

tion' is necessary. The theory of reference 25 does not appear to be 
suitable for use in estimating ^yp^ . 

(2) Fuselage: If the two fuselages are similar in shape, the 

difference in ^Yp^ us can probably be estimated satisfactorily by cor- 
recting for the difference in the relative size of the fuselage and 
wing for the two airplanes. It appears, however, from table X of refer- 
ence 69 unlikely that a reliable prediction of CYpf us can be made 

directly from the geometry of the fuselage. Some additional data on 

Cy A are presented in reference 77- Experimental data from other 
Pfus 

investigations have shown that differences in fuselage cross-section 
can cause very large differences in the variation of ^Yp^ us with 

angle of attack. For example, in the case of a flat fuselage with the 
major cross-sectional axis horizontal, the sign of Cy ft - has been' 

*±us 
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found to reverse at moderate and high angles of attack. Force-test 
data are essential for making estimates in such cases. 

(3) Wing-fuselage interference: For low-wing or high-wing con- 

figurations , wing-fuselage interference causes the value of Cyp to be 

greater than that obtained by adding the contributions of the wing and 
fuselage. (See reference 39-) If the vertical location of the wing on 
the fuselage is generally similar for the two designs, however, any 
correction for a difference in this interference factor can be neglected. 

Vertical tail .- Accurate estimates of ^Yp^ ^ are n - ecessar y 

because this factor is used to estimate the tail contribution to several 
other derivatives. This factor is especially important at low angles 
of attack because in this case the tail contribution is often much 
greater than the wing-fuselage contribution to all derivatives except 
CZp. For this reason it is highly desirable to have tail-off and tail- 
on force-test data for the design under consideration or for a very 
similar design. Corrections to the data for a similar design can be 
made as follows: 

(l) Correction for differences in wing area, tail area, and tail 
lift -curve slope can be made by the following formula: 


( CY Htail) aesign ( Cl Ptall) 


(‘ 


“' a tail 


3 tail 


u ) 


design Sflata 


data ( C I<Hail Stall )data S(iesi e n 


(13) 


The value of CLq.. ... can be obtained from figures 5 an d 6 which are 
^tail 

based on the theory of reference 3^- and on the theory and data of refer- 
ences 28 and 35- The chart of figure 6 can be used to estimate the 
change in the effective aspect ratio of the vertical tail caused by the 
end-plate effect of the horizontal tail. It should be emphasized that 
for the best accuracy the charts in figures 5 and 6 should be used in 
conjunction with formula ( 13 ) Tor correcting existing force-test data 
and not for making a direct estimate of Cyp . . 


(2) In the case of V-tails, the correction for Cy Q can be 

' x Ptail 


made as follows: 


-tail design ( CY|3 V-tail) 


( Kt Vv- 


tail sin 


;„2 r \ 


) 


design D data 


V "tail design \ p V-tail/ da t a /kCt Sy-tail sin : 


J<X N 


sin 2 rj 


data 


5 design 
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where the terms Cl- , r, and K are the same as given in reference 30 
and are defined as follows: 



slope of the tail, lift curve in pitch measured in the plane 
normal to the chord plane of each tail panel 


p dihedral angle of tail surface measured from XY -plane of the 

tail to each tail panel , degrees 

K ratio of sum of lifts obtained by equal and opposite changes 

in angle of attack of two semispans of tail to lifts obtained 
by an equal change in angle of attack for the complete tail 


Values of the term K, which are usually about 0.7, can be obtained from 
reference 30. 


(3) Since large differences in sidewash and dynamic pressure at 
the tail can be caused by differences in wing plan form and wing loca- 
tion, use of experimental data for the specific design or at least for 
a design which has a closely similar wing-fuselage combination and 
vertical tail location is extremely desirable. No methods are available 
which permit accurate predictions of sidewash at the tail, but the 
experimental data of references 39, 49, and 69 can be used to obtain 
some indication of the variation in sidewash with vertical location of 
an unswept wing on a fuselage and the experimental data of references 36 
and 77 provide additional information on sidewash at the tail. Other 
experimental data indicate that the sidewash fields produced by highly- 
swept, low-aspect -ratio wings or by fuselages of flat cross section can 
sometimes be strong enough at high angles of attack to reverse the 
effectiveness of a conventionally-located vertical tail surface. Until 
a reliable method is developed for predicting these large sidewash 
effects, force-test data appear to be the only means by which satisfac- 
tory estimates of ^Yptail can obtained. 


Cnp 

Although attempts have been made to develop methods for estimating 
the yawing moment due to sideslip (static directional stability) deriva- 
tive Cnp (for example, references 68 and 69 ) no reliable method has 

yet been obtained. The use of force-test data therefore seems imperative. 

Force -test data for the design under consideration should be used 
if available. If such data are not available, use data for a similar 
design and correct as explained in the sections to follow. 
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Wing-fuselage . - The corrections for the wing-fuselage contributions 

are: 

(1) Correction for wing - From figure 7 (taken from reference 25) 
the values of (Cn$/ ^L^wing *’ or design under consideration and for 

the design for, which test data are available can.be determined. The 
effect of differences in taper ratio can be neglected. (See refer- 
ences 60 and 66.) The difference between these values of C n ^Ci^ 

should then be added (with proper regard for sign) to the experimental 
data for the complete model. 

(2) Correction for fuselage - The formula 


„ -j ^/Fuselage volume 

%fus ’ ■ 1 - 3 V Sb 

can be used to calculate the C n p of the fuselage (per radian) for the 

design under consideration and for the similar design for which force- 
test data are available. The differences between these two values can 
then be added (with proper regard for sign) to the force -test data for 
the complete model. Formula (15) does not include the effect of fine- 
ness ratio and should not be used for fineness ratios less than 4. This 
formula is an approximate empirical expression which should not be used 

to estimate the value of C nQ directly but should only be used as 

Pfus 

indicated to determine a correction for force-test data. This correc- 
tion method should not be used in the cases of high angles of attack 
when there are large differences in fuselage conf iguration. Force-test 
data are essential in such cases. 

(3) Correction for vertical location of the wing - If the designs 
are generally similar, the correction for the vertical location of the 
wing on the fuselage can be neglected. (See reference 39*) 

(4) Correction for center-of -gravity position - If the center-of - 
gravity position for the design under consideration is appreciably dif- 
ferent from that for the design for which force-test data are available, 
the value of Cnp for the wing-fuselage combination can be corrected 

by multiplying the value of Cyp for the wing-fuselage combination by 

the distance between center-of -gravity positions (expressed in wing 
spans ) . 
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Vertical tail.- Corrections to 


C “Ptail f0r 


C Yp 


tail 


and tail length z/b can be made by the 


differences in 
following formula: 


The 
ences 70 


(c \ < = (c \ ( tail b )desjgn 

\ n ^tail/ design \ ^ tail /data (q y — \ 

V 3tail ^ 'data 

contribution of wing-tip fins to Cnp is treated 
and 84. 


(16) 


in refer- 



In estimates of the rolling moment due to sideslip (effective 
dihedral) Cjp, force-test data for the design under consideration should 

be used. If such data are not available, data for a similar design can 
be used and corrected by the methods that follow. 

Wing-fuselage . - The corrections for wing-fuselage contributions are: 

(l) Correction for wing - From figure 8 (based on reference 25) 
the theoretical values of for the design under consideration 

and for the design for which data are available can be determined. The 
difference between these two theoretical values can then be added (with 
proper regard for sign) to the experimental data. Consideration should , 
be given to scale effect, airfoil section, and surface roughness on the 
value of CZp for highly swept wings. The lift coefficient at which 
the experimental variation of C^ with lift coefficient departs from 

theory is greatest at high Reynolds numbers and for smooth wings with 
round leading edges. For wings with rough surfaces or sharp leading 
edges the effects of Reynolds number on C^g are usually small and low- 

scale wind tunnel data can be used. For airplanes having very smooth 
sweptback wings with rounded leading edges, however, some correction 
should be made for scale effect when estimations are made from low- 
scale wind-tunnel data. Since no rational method has been developed 
for making such corrections it is suggested that, for lift coefficients 
higher than that at which the experimental data departs from the theory, 
an average of the theoretical and low-scale experimental values be used. 
Conservative dynamic stability results will usually be obtained if the 
uncorrected theoretical values of Cjp are used because these values 

are ordinarily greater (more negative) than measured values and because 
the larger negative values of usually tend to decrease the dynamic 

lateral stability. 
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(2) Correction for wing dihedral - The effect of dihedral on Cip 

is treated in references 29, 39, 51, 58, 86, and 79* Correction for 
the difference in dihedral between the two designs can be made by- 
multiplying the incremental geometric dihedral angle (in degrees) by 
the factor obtained from figure 9* A plot of Cjpp against 

aspect ratio for taper ratios of 1.0, 0.5 and 0.25 (obtained from refer- 
ences 58 and 66) and a formula from reference 50 for correcting for 
sweep are presented in the upper portion of figure 9- The lower chart 
and formula in figure 9 (developed from reference 66) should be used in 
addition to the upper chart and formula of figure 9 to estimate the 
values of C 1 ^ for the case of a wing with partial-span dihedral. 

Although this chart and formula apply directly only to wings with one 
dihedral break they can be used to estimate the Cjp for wings with 

two or more dihedral breaks by the method described in reference 66. 

The effect of drooped wing tips and of wing-tip end-plates on Cj R 

p wing 

should be determined by experimental data since no reliable estimation 
procedure for these effects is available. 

(3) Correction for wing -fuselage interference - Although the con- 
tribution of the fuselage alone to C^p is usually negligible, the 

interference between the wing and fuselage can greatly alter the value 
of Cip of the wing. This interference is such that a high location 

of the wing on the fuselage gives more positive effective dihedral 
(higher -CZp) and a low wing location gives less positive dihedral 
than a midwing position. This effect is treated theoretically in 
reference 67 and has been studied experimentally in references 38 to 42, 
The following simplified expression for estimating the increment in Cip 

caused by wing-fuselage interference has been developed from the rela- 
tionships presented in reference 67 and in other sources: 


ac, = 

*3 v b b 


This expression has been found to give reasonably good agreement with 
experimental data for a variety of configurations. It is suggested 
that values of ACjg be calculated from this equation for both the 

design under consideration and for the design for which force-test data 
are available. The difference between these values can then be added 
(with the proper regard for sign) to the force -test data. 


Vertical tail.- The value of 


'3tail 


determined from force -test 


data on a similar design can be corrected as follows to obtain CZg-^-Q 
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for the design, under consideration: 


C2p 


tail 


design 


= Ci 


'hall) 


Ptail ^)design 


data 


z\ 


(18) 


^ Ptail b )<i a -ta 

The results of reference 35 indicate that Ci a can also be affected 

■ Ptail 

by the location of the horizontal tail with respect to the vertical 
tail. If the two designs have approximately the same horizontal tail 
size and location, however, this effect can be neglected. 

The value of ^ 1 3-t a il y ° r a V-ta.il can tie estimated from the 
following empirical formula: 


r- 


( C ^V-tail) 


design \ ^V-tailJ data 


Cy 


^V-tail 


b sin T 


(by. 


tail + W-tail sin *)_ 


Cy PV-tail U ^ 1 j 

b sin T ( b V-tail + Azy-tail sin fj 


design ! 


data 


(19) 


where "by-tail is the developed (not projected) span of the V-tail, 
z V-tail is the vertical distance from the center of gravity to the 
chord of the V-tail (positive up, and T is the dihedral angle of 
the V-tail, More information on V-tails can be found in references 30> 
6l, and 62. 

In the case of a vertical tail located on the wing, there is, in 
addition to the incremental Cjg produced by the tail lateral force, 

an incremental C produced by the interference effect of the vertical 

tail on the wing. Since this interference effect varies greatly with 
spanwise and vertical position of the tail, it should be determined 
from force tests. Usually the interference is such that a vertical tail 
above the wing gives a negative increment of C^p (positive effective 

dihedral) and one below the wing gives a positive increment of C^p. 

In general, the largest interference effects are obtained with vertical 
tails at or near the wing tips. 
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THE ROLLING DERIVATIVES C np , C Zp , C Yp 



The wing and vertical tail are the only airplane components that 
contribute appreciably to the yawing moment due to rolling derivative 
Cnp. The contributions of the fuselage and horizontal tail can usually 

be neglected. 

Wing . - The contribution of the wing to Cnp can be estimated from 

the formula and charts of figure 10 which were taken from reference 86. 
Although these charts apply strictly only to wings having a taper ratio 
of 1.0, experimental data have indicated that they will also provide 
fairly good estimates for taper ratios of 0.50, 0.25 and 0. In the 
estimation formula 



the value of ^Cp 0 j a should be determined, if possible, from force -test 

data obtained at high Reynolds number on the wing under consideration, 
since low Reynolds number data might indicate values of ^CD 0 )a that are 

too large. For the case of smooth wings with a large leading edge 
radius and low or moderate sweep, it is suggested that ^Cp 0 ^ a for the 

airplane be assumed to be zero at all lift coefficients up to the stall. 
This assumption will result in larger negative values of C np than 

would be estimated from low Reynolds number data on ^D 0 ) a con “ 

sequently should lead to conservative dynamic stability results since 
an increase in Cnp in the negative direction has been found to cause 

a reduction in dynamic stability. The value of for highly swept 

wings is often very large at high lift coefficients, especially for 
wings with rough surfaces, sharp leading edges, or triangular plan form. 
For these cases, values of ^D 0 ) a determined even from low Reynolds 

number data might lead to reasonably good estimates of C np . In all 

these cases, however, high-scale drag data should be used whenever it 
.is available. 

Effect of high^lift devices .- The principal effect of leading-edge 
high-lift devices is to extend to a higher lift coefficient the linear 
variation of C np with lift coefficient. The formula and charts of 

figure 10 are directly applicable to this case. The effect of 
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trailing-edge high-lift devices is not so straightforward, hut experi- 
mental data have indicated that the formula and charts of figure 10 
also give reasonably good estimates in this case. 

Vertical tail .- The contribution of an isolated vertical tail 
surface to C n p can be estimated by the following approximate formula 

which has also been commonly used to estimate ^np tail of a complete 

airplane : 


C^Ptail 


- 2 H C *f3 


tail 


( 21 ) 


The values of ^ 

previously discussed. Instead of 
will usually be better to use the 

as determined by force -test data. 


determined from force -test data as 


the geometric tail length 
effective tail length -Cn 


1/ b, it 
^tail / ^tail 


Formula (21) then becomes 


Sail " SStall 


(21a) 


In the case of the conventionally located vertical tail surface, how- 
ever, the rolling wing produces a sidewash at the tail which greatly 
alters the tail contribution to C np . This sidewash causes the values 


of 


Cn Ptail 


to be much more negative than is indicated by formula (21). 


This effect is discussed more fully in reference 36 in which is also 
presented a method for estimating the sidewash. Some preliminary theo- 
retical studies have indicated that the effect of the sidewash on 
C n varies considerably with tail size and tail location and to 


some extent with wing plan form. A comprehensive experimental verifi- 
cation of this theory is planned but as yet only a few scattered checks 
have been obtained. For the case of the conventionally located vertical 
tail surface, the following formula has been found to give estimates 
of Cn that are in fairly good agreement with experimental data: 

Ptail 


or 


Stall 



Stall 


^Ptail 


2 



Cn 0tail 


( 22 ) 


(22a) 
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This formula is based on the assumption that 


Cn Ptall 


is zero at 0° angle 


of attack and varies with angle of attack in the same manner as indicated 

by formula (21). Formula (22) or the method of reference 3 6 can be used 

satisfactorily for first approximations of Cn- for most configura- 

•txail 

tions with conventionally located vertical tails. For more accurate 
estimates, especially for configurations having an unusual tail size or 
tail location, experimental data should be used. 


For wings of triangular plan form with vertical tails either 
directly above or above and slightly behind the wing, experimental data 
have indicated that neither formula (21) nor formula (22) gives an 

accurate estimate of Cn-n but that an average ' of the values obtained 

^tail 

by the two formulas provides a fairly good estimate. 


It is obvious that these methods of estimating Cnp are only 

approximate and are open to question in many cases. Experimental and 
theoretical studies are currently being made to provide better methods 
of estimating Cnp^ a; ^ and, when these methods become available, the 

approximate methods presented herein should be discarded. At the present 
time, however, formula (22) and reference 36 will usually provide much 
more accurate estimates of ^np^ ^ than formula (21) which has been 

in common use up until this time. 



Wing -fuselage . - Most of the rolling moment due to rolling (damping- 
in-roll derivative) C of an airplane is produced by the wing. The 

effect of the fuselage can be neglected unless the ratio of the diameter 
of the fuselage to the wing span is relatively large (greater than 
about 0-3)- For large values of this ratio, the value of C^p will be 

smaller than that for the wing alone by an amount that can be estimated 
from a consideration of the area and lateral center of pressure of the 
wing area included within the fuselage. (See references 103, 108, 
and 112 . ) 

Wing. - The damping in roll of wings has been the subject of many 
experimental and theoretical investigations. (See references on CZp 
in table II.) As a result, some methods of estimating C^p have been 

developed which have been found to give reasonably good agreement with 
experimental results. The method presented in reference 79 appears to 
give sufficiently accurate estimates of C^p for zero lift.- This 
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method is extended in reference 89 to permit the estimation of Cj 

over the normal flight range of lift coefficient. Estimation charts 
and formulas from reference 89 are presented in figure 11. 

High-lift devices .- Experimental data have indicated that the 
damping in roll of wings at low and moderate lift coefficients is not 
greatly affected by the addition of high-lift devices such as trailing- 
edge flaps, leading-edge flaps, slats, and slots. The principal effect 
of such devices is to increase the lift coefficient at which the sharp 
decrease in Ci p occurs. The charts and formulas of figure 11 can be 

used to estimate the C2p of wings with either full-span or partial- 

span high-lift devices with fair accuracy despite the fact that the 
method is not strictly applicable to partial-span high-lift devices, 
(See reference 89 .) 

Wing-tip fuel tanks .- The use of wing-tip fuel tanks usually 
increases the damping in roll of the wing. The experimental data of 
reference 91 for unswept wings indicate that the magnitude of the 
increase varies with angle of attack and depends upon the wing taper 
ratio and on the size and location of the tanks. Unpublished experi- 
mental data indicate similar effects of wing-tip tanks on sweptback 
wings. The following approximate formula for estimating the increment 
in C^p produced by wing-tip tanks at low lift coefficients is based 

on the limited amount of available experimental data and should not be 
expected to yield very close quantitative estimates: 


AC r 


£7 tanks 


= Cl 


tanks off 




Maximum tank diameter \/ v / ^ 

T.T-i m rr o-no m J ' ^ ' \ J 


Wing span 


where, for symmetrically mounted tip tanks, 

k t = 6 

for tanks mounted below the wing tip or forward on the wing tip, 

K T = 3 

and for pylon-mounted tip tanks, 

= 1 

Experimental data for both unswept and swept wings indicate that 
f AC7pj^ an ^. s usually becomes smaller with increasing angle of attack 

and, in some cases, actually reverses sign at high angles of attack so 
that the tanks are decreasing rather than increasing the damping in roll. 
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The data of reference 91 can he used to obtain an approximate estimate 
of the effect of angle of attack for unswept wings. 

Tail surfaces . - The contribution to Cjp of conventional type 

horizontal, and vertical tail surfaces is usually very small and, in 
most cases, negligible. When an airplane rolls, the wing produces a 
rotation of flow at the tail surfaces which reduces the already small 
damping moments of the isolated surfaces, except in the case of the 
vertical tail at high angles of attack where the tail center of pressure 
is below the center of gravity. 

The contribution of an extremely large horizontal tail to C 2 p 

might not be negligible and can be estimated by multiplying the value 

of C 7 for the particular tail plan form obtained from the charts and 
p S t 7b t \2 

formulas of figure 11 by the factor 0.5 in which the fac- 

tor 0.5 is included to account for the rotation of flow produced by the 
wing . 

i 

The contribution of an isolated vertical tail surface to is 

given by the following approximate formula: 

C! Ptall = 2 (sf ^tail (2k> 


As in the case of this formula can be modified to provide an 

approximate correction for the effect of the wing on the damping in 
■roll of conventionally located vertical tail surfaces: 


-(IU 

An analysis of this expression indicates that 

negligible at .low and moderate angles of atta< 
but that it might be fairly important at very high angles of attack 
where z/b is a large negative value. As in the case of Cn p , experi- 
mental data indicate that, for a vertical tail located either directly 
above or slightly behind a wing of triangular plan form, the value of 
C Zpt ^ can be .estimated with better accuracy by an average of 

formulas (24) and (25) than by formula (25) alone. For conventional 
tail arrangements, however, formula 25 gives better correlation with 
experimental data. 


CYp 


tail 


(25) 


the value of C 7 is 

Ptail 

:k where z/b is positive 
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Wing . - The following formula for the derivative Cy-p (lateral 

force due to rolling) from reference 86 is based on experimental data 
and is the same as that presented in reference 25 except for an addi- 
tional correction to account for tip suction: 


cl 


A + cos A 
A + 4 cos A 


tan A + — 
A 


(26) 


The data of reference 86 show that this formula applies only for lift 

cl 2 

coefficients below that at which the drag factor Cp - begins to 

increase. At higher lift coefficients the experimental data indicate 
smaller values of Cyp than given by formula (26). For these cases an 

approximation of the value of Cyp can be obtained from the experi- 
mental data of reference 86. As in the case of Cn p , the break in the 
variation of Cyp with lift coefficient should be expected to occur 

at lower lift coefficients for wings having sharp leading edges or 
rough surfaces and for wings tested at low Reynolds numbers. 


Vertical tail.- The discussion concerning Cn^ and CZ-n 

Ptail Ptail 

is also applicable to Cy-n * The value of Cy_ for an isolated 

-Ptail -Ptail 

tail surface is given by the formula; 


C Y P = 2te)cy 3 • (27) 

p tail \ 7 p tail 

This formula can be modified as follows to account approximately for the 
effects of wing sidewash in the case of a conventionally located vertical 
tail : 



tail 





tail 


(28) 


An average of formulas (27.) and (28) can be used for tails located 
either directly above or above and slightly behind the wing. 
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THE YAWING DERIVATIVES C nr , Cz r , AND Cy r 

Cn r 


Wing-fuselage . - In the past, the contribution of the wing-fuselage 
combination to yawing moment due to yawing (damping in yaw) derivative 
Cn r has usually been found to be small compared to the contribution of 

the vertical tail. The fuselage contribution to the damping in yaw 
depends, of course, on the relative size of the fuselage and wing. In 
the past, the relative size of these components has generally been such 
that the fuselage contribution could be neglected. (See references 82 
and 83.) For some recent designs which have a large I'uselage relative 
to the wing, however, the fuselage contribution to C nr is important. 

In the case of fuselages having flat sides or having a flattened cross 
section with the major axis vertical the fuselage contribution may also 
be important and some fuselage contribution to Cn r should be assumed, 

especially at high angles of attack. On the other hand, experimental 
data have shown that a flattened cross-section fuselage with the major 
axis horizontal can have negative damping in yaw at moderate and high 
angles of attack. 

The contribution of the wing to Cnr can b e estimated from the 
formula and charts of figure 12 which were taken from reference 25 . 
Values of Cp 0 for the wing should be estimated from force-test data. 

For values of x/cf greatly different from zero, the charts of refer- 
ence 25 can be used. The formula and charts of figure 12 are not con- 
sidered reliable at high angles of attack, especially for swept wings. 
The use of experimental data from the references on Cn r listed in 

table II is recommended in this case. 

The effect of partial-span inboard flaps on C nr can usually be 

neglected. (See reference 82.) The effect of full-span trailing-edge 
or leading-edge high-lift devices can be estimated satisfactorily from 
the formula and charts of figure 12. Values of Cd 0 in this case are, 

of course, for the wing with the high-lift device installed. 

Vertical tail .- The contribution of a conventional -type vertical' 
tail to Cn r can be estimated from the formula 


C "rtaU ' 2 (5f C *P tall 


(29) 
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or, with the effective tail length 
the geometric tail length i/b. 


■ Cn Hau/ CY Stall substltuted f ° r 


'nr. 


tail 


J Cn Ptaix) 2 

' ^tail 


(29a) 


The experimental values for C nT . presented in reference 82 for 

Hail 

power-on or flap-down configurations are 30 to bO percent greater than 
values predicted by formulas (29) or (29a). These differences are 
attributed to lag of sidewash effects in the free-oscillation tests 

used in measuring Cn r . In estimations of C nT , for stability cal- 

' 1 I tail 

culations, similar lag of sidewash effects should be assumed if the 
oscillatory mode is of primary importance but no lag of sidewash should 
be assumed if the aperiodic mode is most important. 


Methods for estimating the c n r ^ ail 
are presented in references 70 and 82. 


for wing -tip vertical tails 


Cz 


r 


The wing and vertical tail are the only airplane components that 
contribute appreciably to rolling-moment -due -to-yawing derivative Cz r 

of an airplane. The contributions of the fuselage and horizontal tail 
can usually be neglected. A semiempirical method for estimating Cz r 

is presented in reference 85. This method involves the use of experi- 
mental data on the parameter Czp to correct the theoretical values of 

Cz r wing § iven reference 25 and to estimate the value of Curving* 


Wing . - The formula of reference 85 and the charts of Cz^Cl from. 

reference 25 for estimating Cj- ■. are given in figure 13. The 

x wing 

values of Cl to be used in the charts can be obtained from fig- 
ure 8. For taper ratios less than 0.25, values of and Oz^/Cl 

for a taper ratio of 0.25 can be used. The value of Czr used in 

the formula should be the same as the value of C 


3exp 
estimated from 


^wing 


experimental data by the method indicated in the section on Cz^. In 
the case of Cz r , however, (unlike the case of CZp) conservative 
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dynamic stability results will usually be obtained if the smaller values 
of the derivative (based on low-scale experimental data) are used 
instead of the larger (theoretical) values. This difference is a result 
of the fact that either an increase in the normally negative value of 
Cjp or a decrease in the normally positive value of Cz r can cause 

reduction in dynamic stability. As pointed out in reference 85 the 
estimation procedure shown in figure 13 appears to account satisfactorily 
for the effects of high-lift devices, wing, dihedral, and airfoil section. 

Vertical tail .- The contribution of the vertical tail to C z r is 
usually estimated by the formula 

C *rtail ' -KsXffcutall (3 °> 


where CY Ptall 


is preferably obtained from force-test data. 


When 


experimental data on *C are available, the following formula from 

reference 85 can be used and will probably be more reliable than equa- 
tion (30) because it takes into account any interference effects that 
might cause the effective vertical location of the center of pressure of 
the tail to be different from the location determined by geometrical 
procedures: 


Cz 


r tail 


= -2 f C 


b/ ^Ptail 


(31) 


or with the effective tail length 
geometric tail length Z/b, 


^ n 0tail/^3tail 


substituted for the 


f c n, 


Cz = 2 

r tail 


Ptail' 


cy 3 


Cz 


tail' 


p tail 


(31a) 


C *r 

Wing . - The theory of reference 25 gives values of the derivative 
Cy r (lateral force due to yawing) for the wing for a taper ratio of 1.0. 

The experimental data of references 25 and 59 indicate that this theory 

is inadequate for making reliable estimates of Cy~ . It is recom- 

r wing 

mended therefore that the experimental data given in references 25, 58, 
59, and 60.be used in making estimates of ^Yr W i n • 
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Vertical tail.- The value of 'Cy,. can be estimated by the 

■ . ^tail 


formula 


CYr tail ‘ 2 b ° Y Ptail 


( 32 ) 


or by the formula in which the effective tail length 
is substituted for the geometric tail length z/b: 


" Cn Ptail/ CYp tail 


CYr tail 2Cn Ptail 


(32a) 


The discussion of lag-of-sidewash effects for 
C ^ r tail* 


Cn-r apply also to 

r tail 


EFFECTS OF MACH NUMBER 


The effects of Mach number on the lateral stability derivatives 
have been treated theoretically in many investigations (see table II) 
but very little experimental data have been obtained to verify this 
theoretical work. Moreover , only a small part of this experimental 
work has been covered in published reports (reference 111) because most 
of it is classified at the present time. It appears, therefore, that 
estimates of the lateral-stability derivatives for the time being will 
have to be based largely on theoretical work. 

The effects of Mach number on the stability derivatives can be 
usually considered negligible for all airplane components except the 
wing and vertical tail. For the low-lift-coefficient condition in the 
case of many high-speed airplanes, the vertical tail contributes more 
than the wing to all the stability derivatives except CZp. For this 

reason, in calculations for transonic or supersonic speed conditions it 
is especially important to know the effects of Mach number on the 
vertical-tail lift-curve slope or Cy^ ^ 1 . 

Wing . - The effects of compressibility on the subsonic stability 
derivatives of the wing can be estimated by the formulas of reference 26. 
The values of the supersonic stability derivatives for some wing plan 
forms can be estimated by the references tabulated in table II. In this 
table the derivatives are grouped according to the type of wing plan 
form and to the particular derivatives covered. A helpful summary and 
discussion of the effects of Mach number on the derivatives for several 
different wing plan forms is presented in reference 103 • A summary of 
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the theoretical lift -curve slope , damping in roll, and center-of- 
pressure characteristics of various wing plan forms is presented in 
reference 107 . In the cases in which the theory shows large or abrupt 
changes in a stability derivative with changes in Mach number (for 
example, fig. 10 of reference 103 ) special care should be taken in 
estimating the derivative in that particular Mach number range. The 
abrupt changes should be smoothed or faired out in a manner similar to 

that suggested in the following section for estimating Cy n 

p tail 

In some cases, experimental data for supersonic speeds will be 
available on the sideslip derivatives and on the damping-in-roll deriva- 
tive C^p. In such cases the experimental data should be used in pref- 
erence to the theory. Some experimental results have indicated that the 
effect of the vertical location of the wing on the fuselage on the 
derivative C^p might be greatly different at supersonic speeds from 

that at subsonic speeds. Since no methods are presently available for 
estimating this effect for the supersonic case, it appears that, at 
least in the case of high -wing and low -wing designs, force -test data 
are necessary for obtaining an accurate estimate of C^p. 

Vertical tail .- The sideslip derivatives produced by the vertical 
tail at transonic and supersonic speeds can be estimated theoretically 
but should be obtained from force-test data whenever possible. These 
sideslip derivatives can be used to estimate the tail contributions to 
the other derivatives as pointed out previously. In estimates of the 
value of Cyp^ ^ for transonic and supersonic speeds, corrections must 

be made for the effect of Mach number on the lift-curve slope of the 
tail, and these corrections should account for any differences in the 
end-plate effect of the horizontal tail on the vertical tail. 

For Mach numbers below about 0.8 or 0.9 and above about 1.6 or 1.8 
the effect of Mach number on the lift-curve slope of the vertical tail 
can be estimated satisfactorily from the theoretical values of refer- 
ences 2 6 , 34, and 107. Since experimental data indicate that theoreti- 
cal values of lift-curve slope are usually too high for Mach numbers 
from about 0.8 or 0.9 to about 1.6 or 1.8, the empirically determined 
fairings shown in figure ih are recommended for use as a guide in the 
use of the theory to obtain approximate estimates in this Mach number 
range when force -test data are not avilable. 

Experimental data have indicated that for vertical-tail configura- 
tions which have a tail length (distance from the center of gravity to 
the tail center of pressure) that is relatively short in terms of tail 
chords, the rearward shift of the tail center of pressure at supersonic 
speeds can cause an appreciable increase in the tail length and 
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consequently an appreciable, increase in the magnitude of some of the 
tail derivatives. Theoretical center-of -pressure positions for various 
plan forms, at supersonic speeds are given in reference 107 • 


EFFECTS OF POWER 


On the basis of existing information, the effects of power on the 
lateral stability derivatives appear to be negligible in the case of 
jet-propelled airplanes but these effects are often very large in the 
case of single-engine propeller-driven airplanes. Methods are available 
for estimating some of these power effects but in most cases experimental 
data are necessary for making a satisfactory estimate. The effects of 
power can be broken down into two general classes: 

(1) The effects of the lateral fprce produced by the propeller 
itself 


(2) The effects of the propeller slipstream on the wing, fuselage, 
and vertical tail of the airplane 


Effects of propeller lateral force . - A method of estimating the 
propeller-lateral-force derivative Cy^ is presented in reference 31 

which is based on the work of references 32 and 33* The contribution 
of the propeller lateral force to the other stability derivatives can 
be estimated from this derivative by assuming that the propeller is 
effectively a vertical tail surface and by using the expressions for 
the tail contribution to the various derivatives presented in the 
preceding sections. Some experimental data on the effect of windmilling 
propeller on all of the derivatives are presented in reference 65 . 

Effects of propeller slipstream .- The effects of propeller slip- 
stream on the lateral- stability derivatives are usually much greater 
than the effects of propeller lateral force in the case of single-engine 
tractor airplanes. The slipstream effects on the wing, the fuselage, 
and the vertical tail can be considered as three independent effects. 


The slipstream effects on the 
for the derivatives C 2 p and Ci r 

decrease in effective dihedral ^-C 

planes are presented in references 
highly desirable to determine this 
because interference effects make 
very difficult. The effect of the 


wing can usually be neglected except 
.v Experimental data showing the 

with power for single-engine air- 

55; 56, 7^-; and 80. It appears 
effect of power experimentally 
accurate estimations of the effect 
slipstream on the value' of Cj. 


cannot be estimated from the data on 


Pwing 


tr wing 
as described in the 
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section on C i T . In fact, this procedure would probably give the wrong 
sign for the increment of c ^r v i n g contributed by the slipstream. An 

approximation of this increment might be obtained by estimating the slip- 
stream velocity and the lateral displacement of the slipstream caused 

by yawing. Usually the power effects on arid will be 

p wing r wing 

greatest for the flap-extended configuration. 

In the case of the single -engine airplane the effect of the slip- 
stream on the fuselage is usually to increase negatively the values 
of Cnp and Cyp. (See references 5^, 55, 5 6 , 71* 7^, and 7 6 .) Since 

no accurate methods of estimating these slipstream effects on Cnp 
and Cyp are available, it is necessary to determine them from force - 
test data. 


The effects of the slipstream on the vertical tail are often very 
important and should also be determined from experimental data, if 
possible. The- increase in dynamic pressure at the tail caused by the 
slipstream is treated theoretically in reference 116 and is illustrated 
by the experimental data of references 50 , 5 ^, 55 , 56 , 71 , 7 ^, and 76 . 
The experimental data of reference 76 also show that the propeller slip- 
stream can cause a destabilizing sidewasji at the tail which will tend 
to reduce the stabilizing effect of the increased dynamic pressure at 
the tail. Since these data indicate that slipstream effects on the 
vertical tail vary greatly with airplane configuration and propeller 
arrangement (single or dual rotation), use of experimental data appears 
to be the only satisfactory estimation procedure at present. 

Suggested estimation procedure for power effects .- The following 
procedure is suggested for estimating power effects. Obtain force- 
test data for tail off and tail on. Use tail-on data directly for Cyp, 

Cnp, and C^p. Estimate rolling and yawing derivatives as follows: 

(l) Estimate ^^ ro ^ e ^± er ^ rom reference 31 and use this deriva- 
tive and proper linear dimensions to estimate other propeller deriva- 
tives (rolling and yawing derivatives) in the same manner as tail 
derivatives . 


c Yp 


(2) Subtract tail -on data from tail -off data to get values of 


tail' 


Cn Ptail' and 


Cl 


Ptail 


for the power-on condition and use 


these values to estimate the tail contribution to the other derivatives. 


(3) For tail-off values of rolling and yawing derivatives, use 
same values as for power-off for all derivatives except C i T . Estimate 

C i r as suggested in preceding section. 
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(4) Add the values obtained in steps 1, 2, and 3 to get the 
rolling and yawing derivatives for the complete airplane. 


INADEQUACIES IN PRESENT INFORMATION AND METHODS 


In the course of summarizing the estimation methods for the various 
stability derivatives, the need for much additional information on all 
the derivatives became apparent. In particular, information is needed 
to aid in the estimation of the derivatives in the transonic and super- 
sonic speed ranges. Additional work also needs to be done in correlating 
and analyzing existing subsonic data and in obtaining new experimental 
data for the development of semiempirical methods of estimating the sub- 
sonic derivatives without resort to force-test data. Another important 
need is 'for full-scale experimental results at all speeds for checking 
both low-scale data and the existing methods of estimating derivatives. 
Details of the need for additional work along these lines are discussed 
in the following sections. Studies should also be made to determine 
the conditions for which the use of steady-state stability derivatives 
in conventional stability equations is inadequate and to determine 
satisfactory methods of treating such conditions. 


Transonic and Supersonic Speeds 

Additional theoretical work is needed on the estimation of sta- 
bility derivatives in the transonic and supersonic speed ranges to 
cover the range of wing plan forms for all the derivatives. In particu- 
lar, more work is needed on plan forms currently under consideration, 
such as wings having moderate sweepback and taper. This need is illus- 
trated by table II which indicates that very little material is available 
on the stability derivatives for such plan forms except, perhaps, for 
the derivative C^p. It appears from the table that this derivative 

and the triangular plan form have, in the past, received a dispropor- 
tionate .share of attention, probably because of the greater ease with 
which they could be treated theoretically. 

The greatest need for work on stability derivatives at the present 
time is probably in the measurement of the derivatives at transonic and 
supersonic speeds. Experimental data on wings are urgently needed for 
checking the theoretical work and for use in the development of empirical 
corrections to the theory wherever necessary. Such corrections are 
particularly needed for fairing out abrupt variations of the derivatives 
with Mach number and for fairing through the Mach number range for which 
theory predicts infinite values. Examples of such discontinuities as 
indicated by theory are shown in figures 8 to 13 of reference 103 • 
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Since experimental data obtained at supersonic speeds on ving-fuselage 
combinations and on complete models have revealed interference effects 
that are different from those obtained at subsonic speeds, it appears 
highly desirable to obtain at least a- limited amount of experimental 
data at transonic and supersonic speeds to evaluate these interference 
effects. For example, investigations should be undertaken to determine 
the effect of wing-fuselage interference on the derivative C and 

the end-plate effect of the horizontal tail on the lift -curve slope of 
the vertical tail. 

Most of the experimental data on stability derivatives at transonic 
and supersonic speeds will of necessity be obtained at Reynolds numbers 
considerably less than full-scale values and under test conditions which 
might render the results open to question in some cases. Full-scale 
checks in flight of the low-scale data and of the estimation methods 
therefore appear to be desirable. Consequently the methods of measuring 
stability derivatives in flight now being developed by the Cornell 
Aeronautical Laboratory, the Massachusetts Institute of Technology, and 
the NACA should be extended to transonic and supersonic speeds when the 
methods appear to be developed to a satisfactory degree of reliability 
for the subsonic case. Some preliminary considerations involved in the 
use of these flight techniques are discussed in references 117 to 120. 


Subsonic Speeds 

The methods presented in this paper for estimating the stability 
derivatives at subsonic speeds depend either directly or indirectly on 
the use of force -test data. These methods are probably more reliable 
than methods which do not involve the use of force -test data on the. 
particular design under consideration or on a similar design. Methods 
which do not rely on such data are desirable in some cases, however, 
because the necessary data will not always be available. 

In the case of sideslip derivatives, empirical methods can probably 
be developed largely from existing information. In some cases it will 
be necessary to augment the existing information with new results since 
much of the available force -test data were not obtained in a manner that 
would make the data readily usable for developing general estimation 
procedures. 

In the case of rolling and yawing derivatives, considerably less 
information is available than in the case of the sideslip derivatives. 
Most of the information now available was obtained in the Langley 
stability tunnel, principally on wing conf igurations and to a limited 
extent on complete airplane models and airplane components other than 
the wing. Considerably more work is required, especially for components 
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in combination, before satisfactory methods can be developed for 
estimating rolling and yawing derivatives without the use of force-test 
data on the particular design under consideration or on a similar design. 

In discussing the work necessary for developing new procedures for 
estimating the stability derivatives without the use of force-test data 
on the design under consideration or on a similar design, it is useful 
to break the problem down into two parts: (1) effect of individual 

components and (2) the effect of interference of the components on each 
other. 

The principal components to be considered are the fuselage, wing, 
vertical tail, and propeller. For the isolated fuselage, the main 
problem is the development of methods for the estimation of C n p and 

then, perhaps, of Cn r and Cyp. For the isolated wing, the main 

problem is to estimate the derivatives at lift coefficients above that 
at which separation begins. Such estimations can be made with reasonable 
accuracy for some of the derivatives- by existing methods which make use 
of force-test data, but the development of methods which do not involve 
the use of force-test data will probably be very difficult. For the 
isolated vertical tail, the problem is to establish the effective tail 
area and aspect ratio from the geometry of the tail so that the lift" 
curve slope (or Cy^ of the tail can be calculated. Solutions to this 

seemingly simple problem have in the past become involved with interfer- 
ence effects so that, as yet, no reliable methods have been published 
for estimating Cyp of the vertical tail from its geometry. For the 

isolated propellers, the work that is needed at present is a systematic 
check of existing methods of estimating the lateral force on the 
propeller to determine the accuracy of these methods. 

The principal interference effects to be considered are mutual 
interference of the wing and fuselage; wing-fuselage interference on 
the vertical tail; horizontal -tail interference on the vertical tail; 
propeller-slipstream interference on the wing, fuselage, and vertical 
tail. The mutual -interference effects of the wing and fuselage are 
probably important only for the derivatives Cjg, C n p> an< ^ Cj r ' ^ 

large amount of experimental data is available for the sideslip deriva- 
tives but no procedures for estimating the interference effects on these 
derivatives have been reported. Wing-fuselage interference has very 
important effects on Cyp of the vertical tail, and consequently on all 

of the stability derivatives for some flight conditions. These effects 
result from the sidewash and change in dynamic pressure at the tail 
which may result from sideslipping, rolling, or yawing. Although con- 
siderable data which show these interference effects are available, 
particularly for the case of sideslipping, no reliable methods exist 
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for estimating the interference effects. Horizontal -tail interference 
also has an important effect on Cyp of the vertical tail for some 

horizontal -tail positions. Some vork on a limited number of configura- 
tions has been done toward developing methods of estimating this effect 
but data are required on more configurations before the generally appli- 
cable methods can be evolved. The propeller slipstream can cause impor- 
tant effects on Cjp and C i T of the wing, on C n p and Cyp of the 

fuselage, and on Cyp of the tail (and consequently on the tail contri- 
bution to all the derivatives). Some data are available for the effect 
of the slipstream on the sideslip derivatives but, because of the com- 
plexity of this problem, considerable additional data may be required 
before a satisfactory method of estimating the slipstream effects- can 
be developed. 

As mentioned in the preceding section, full-scale checks of low- 
scale data and of the estimation methods are desirable. For the sub- 
sonic case some of the checks can be obtained from large-scale wind- 
tunnel tests but some checks in full-scale flight tests should also be 
obtained when the various methods of measuring stability derivatives in 
fli u nt have been developed to a satisfactory degree of accuracy. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va., December 13* 1950 
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APPENDIX A 


EQUATIONS OF MOTION 


The dimensional equations for the lateral motions of an airplane 


are 


.2 A dL d0 , d% dL dt _ 

mkv ^ + ^2 " It - v - L c = 0 


"x dt 2 - dt 


(Al) 


, d H dN d0 ^ , 2 <^t dN dt dN _ 

^xz H| - ^ dt + ^z 2 ^2 - SF dt " & v - N c " 0 


( A2 ) 


"Ip It " ( Lift )0 + mV H " S f+T - (Lift) (tan 7)* + m §J 


dt 5r dt 


<3Y 

dt ~ 


v - Y c = 0 


(A3) 


If equations (Al) and (A2) are divided .by ipV%b and equation (A3) is 

1 O 

divided by s th e equations of motion may be expressed in the con- 

ventional nondimensional form in which they have generally been presented 
in NACA reports (for example, see reference 2): 


W 0 - | %-g * 2^ xz 0 - i c lr f - - C !c = 0 


fi.ic M 


2 d__t 1 p dt 
ds 2 


^%Z - I c n p di + 2liK Z ' 5 c »r S - V - <V - ° 


> (Alt) 


I % if - + 


<V - Cy c = 0 


dt 

ds 


I ° Y r is " C L( tan 7)t + 2|i i| 
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In order to convert these equations into a form which will reduce the 
number of arithmetical and algebraic steps in performing stability cal- 
culations, equations (A4) are multiplied by m/pSb and written in the 
following form: 


D 2 - ZpD)0 + ^K]D 2 - ZjJ))* - ZpP - Z c = 0 
(KgD 2 - UpE)) 0 + ^D 2 - n^Dj* - n p p - n c = 0 

(' y P D " + (d - y r D - ^ tan y)i|r + (d - y p )p - 


^ (A5) 


y c = 0 


where 


H = 


m 

pSb 


T ~ 


m 

pSV 


t 

a = - 


D - A 

da 


K xz 

1 ■ 72 

K x 


K 


K 2 = ,, 2 


XZ 


= ~^~2 C i B 
2K V P 


n R = 


2K< 


2 C n p 


y 3 " 2 C Y F 


l v = C Z 

p 4 k v 2 P 


Hp 


** 8 "p 


y P = 1ST c x 


P 


h = p C h 

4K y 


n r = 


4K, 


2 n r 


Yr = 




lr = 


2K-< 


2 ° z c 


n r = 


= -JL 
2K- 


2 ~n c 


7c - 2 Cy c 
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' APPENDIX B 

APPLICATION OF THE LAPLACE TRANSFORMATION TO CALCULATING MOTIONS 


The application of the Laplace transformation to the calculation 
of the lateral motions of airplanes is presented in order to illustrate 
the development of the equations of motion in the form in which they 
are presented in the present paper. This work is similar to that pre- 
sented in references 5 and. 6. In fact, it follows the presentation in 
reference 5 very closely. Reference 6 presents a brief explanation of 
the Laplace transformation and its application to solution of the equa- 
tions of motion of an airplane. This paper also makes reference to 
detailed explanations of the Laplace transformation. In cases where 
modification of the equations presented in the present paper are neces- 
sary, reference should be made to these texts for an understanding of 
the mathematics involved. Applying the Laplace transforms 

L(l) = i L(D0) = X0 X - 0 O 

L (0 ) = 0x L(d20) = X% - \0 O - (D0) q 


and. multiplying each of the equations by X transforms equations (A5) 
from appendix A to 



where 

r x = (X 2 - Z p M0o + (l^X 2 - Z r X)i 0 + X.(D0) o + K 1 \(D^) 0 + l Q 
r 2 = (K 2 k 2 - n p X.)0o + ( X 2 - n r X)^ 0 + K 2 X(D0) O + X(Di) 0 + n c 


r 3 = -y p x 0o + ( x - y r ^o + x Po - yc 
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Solving equations (Bl) by determinants gives 


-ZpX 

r l 

K]X 3 - Z r X 2 

-n p X 

r 2 

X 3 - r^X 2 

4 2 - ypX 

r 3 

X 2 - y r X 2 - ^(tan y)X 

-ZpX 

x 3 - z p x 2 

K l X 3 - z r x 2 

-n n X 

K 0 X 3 - n X 2 

X 3 - n X 2 

(3 

2 P 

r 

X 2 - y p X 

,2 C B 
" y P X ~ ~ X 

X 2 - y r x 2 — J^tan /)X 


which may be expressed as 


0) 


agX^ + a^X^ + agX 3 + a^X 2 + a^A. + a^ 
X 2 (ax 4 + BX 3 + CX 2 + DX + e) 


(B2) 


Similarly, the expressions for \|r^ and |3^ are 


^X 


b 0 x5 + bj_X^ + bgX 3 m- b^X 2 + b^X + b^ 


! ( 


AX^ + BX 3 + CX 2 + DX + E 


I 


(B3) 


*X 


CqX^ + C-]_X 3 + CgX 2 + C^x + 

X ^AX^ + BX 3 + CX 2 + DX + e) 


(B4) 


where the expressions for the coefficients in equations (B2) to (B4) are 
given in terms of the mass and aerodynamic stability derivatives by 
equations (1) to (4) in the main body of this paper. 


In order to obtain the actual variables from the transformed 
variables, an inverse Laplace transformation must be applied. The 
expressions for 0^, and |3^ are of the form where u^ 
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and v^ axe polynomials , the degree of being higher than that 

of u^. The inverse transform of a function of this type is 


T-i/M = e^n 

\v\) 5=1 v'(x n ) 


(B5) 


In this equation all of the roots X of = 0 are assumed to be 
distinct. This assumption is valid for 3\J but for 0^ and \|r^, 

v^ = 0 has two zero roots. (See equations (B2), (B3), and (B4).) The 
terms in the equations for 0 and \|r resulting from the two zero roots 
are 


where 


§(0) + fi(0)o 


a = — x 2 
v x 


(b6) 


The inverse transforms of 0^, ^ and 3\ are from equations (B5) 

and (B6) 

0 = A^e°^l + A2e a ^2 + A^e 0 ^ + A^e 0 ^4 + + Ag (B7) 

4 

= B-^e^l + B2e a ^ + B^e 0 ^ + B^e 0 ^ + B^a + Bg (B8) 


3 = C^e^^ + + c^e°^ + 


(B9) 


The equations for the rolling velocity p and the yawing velocity r 
can be obtained from equations (B7) and (B8) by differentiation 

p .= Tjr(-AqX-^e ^ + AgXge 2 + A^X^e 3 + A^Xj^e 4 + A^ (BIO) 




ctX-i 


OXp 

+ BgXge + B^X^e 


oX- 


3 3 


B^X^e 


aX^ 


+ B 


(Bll) 


where the expressions for the coefficients of equations (B7) to (Bll) 
are given by equations (6) to (8) in the section "Calculation of Motions." 
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APPENDIX C 

SOLUTION OF BIQUADRATIC EQUATION 


Many methods are available, of course, for solving for the roots 
of a biquadratic equation. For example, there are Horner's, Ferrari 's, 
Bernoulli's, Descartes', and Hitchcock's methods; various methods of 
solution by trial; and also various graphical methods such as that 
given in reference 1. Solution by trial in which synthetic division is 
used, however, is recommended as being the simplest method for most 
lateral stability work. The characteristic equation for the lateral 
motions of an airplane 

AX^ + B\ 3 + CX 2 + D\ + E = 0 

generally has two real roots and a pair of conjugate complex roots. For 
these cases the two real roots can be factored out easily and the 
remaining quadratic solved for the conjugate complex roots. In the few 
cases for which all four of the roots of the characteristic equation 
are complex, Descartes* method can be used to factor the biquadratic 
equation into two quadratics. When there are real roots, solution by 
Descartes* method requires more time than factoring out the real roots 
singly and consequently is not recommended for general use. These 
methods of solution are explained in the following sections. 


Solution by Trial by Means of Synthetic Division 

Solution for real roots by trial by means of synthetic division 
consists of successive approximations of a root and checking by syntheti 
division until the root is determined to the desired degree of accuracy. 
This check by synthetic division is based on the fact that if a is a 
root of a polynomial f(x) then x - a is a factor of f(x) and con- 
sequently no remainder is left when f(x) .is divided by x - a. 

The method of solving the stability biquadratic equation by trial 
with synthetic division is explained in three steps in the following 
sections. First, the rule for synthetic division and a numerical 
example are given. Second, the specific use of synthetic division for 
factoring a biquadratic is illustrated by a simplified example for 
which the roots are known. This example shows how the cubic and then 
the quadratic factors of the biquadratic are obtained. Third, the use 
of synthetic division in extracting the roots of a representative 
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characteristic stability biquadratic is illustrated with special refer- 
ence to methods of making the first approximations of the real roots. 


Explanation of synthetic division .- Synthetic division is explained 
in almost all algebra text books but is presented herein for the con- 
venience of the reader. The rule for synthetic division may be given 
as follows: 

Assume that a polynomial in x (f(x)) is to be divided by x - aj 
write the coefficients of the polynomial in order, supplying 0 when 
a coefficient is lacking. 

Multiply a by the first coefficient, and add (algebraically) the 
product to the next coefficient. 


Multiply this sum by a, add to the next coefficient, and proceed 
until all the coefficients are used. The last sum is the remainder and 
also the value of the polynomial when a is substituted for the 
variable x. 

4 3 2 ^ 

For example, divide x + 3 X + 3 X - x- 6 by x - 3 - 


1 + 3 + 3 - 1 -. 6 


+ 3 + 18 + 63 + 186 
1+6+21+62+180 


3 


Use of synthetic division in factoring out roots . - The use of 
synthetic division to factor out two known rational roots of a biquadratic 
equation is illustrated by the following simple example. These two 
rational roots represent the two real roots of the characteristic sta- 
bility equation which, of course, are not normally known but can be 
approximated by the method given in the next section of this paper. 


One factor of the biquadratic is x - 1 so there is no remainder 
when the biquadratic is divided by the root 1 


1+3+3-1-6 


+ i + ij- + 7 + 6 

1 + 4 + 7 + 6 0 


1 


Since the remainder is 0 , x - 1 is one factor of the biquadratic 

equation and x^ + 4 x^ + 7 X + 6 is another factor. Inasmuch as a 
cubic equation must have at least one real root, a second real root of 
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the biquadratic equation can be factored out of the cubic . For 
example x + 2 is a factor so divide the cubic by the root -2. 


1 + 4 + 7 + 6 


- 2 - 4 - 6 

1 + 2 + 3 0 


-2 


2 

The factors of the biquadratic then are x - 1, x + 2, and x + 2x + 3 
The quadratic factor can be solved for its roots by the quadratic 
formula. For example 


x 


-2 ± \[k - 12 
2 — 


-1 ± i \[F 


Example of application to characteristic equation .- Reasonably 
accurate first approximations to the real roots of the characteristic 
equation can be obtained from simple formulas. Successively closer 
approximations can then be obtained by interpolating from the remainders 
The following example illustrates the application of this method to 
obtaining the roots of the stability biquadratic. The biquadratic 

X 1 *- + 10.43x3 + 16. 32X 2 + 68. 6x - 9.10 = 0 


is of the form 


^ AX^ + BX 3 + CX 2 + DX + E = 0 

Since the coefficient E is generally much smaller than coefficient D 
in lateral stability work, one of the real roots (usually the smaller 
of the two) is approximately equal to -E/D or it may be more closely 
approximated by the equation 




D 


CE 

D 


or for the particular case 


-9.10 

7TZ ( 16. 32 ) ( -9.10) 

60 - 6 60 


0.129 


x 
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Approximating the root by synthetic division 


1 + 10.43 + 

16.32 + 

68.6 - 9.10 


Approximation 

+ .13 + 

1.36 + 

2.3 + 9-10 

.1284 

2 

+ .13 + 

1.36 + 

2.3 + 9-14 

.129 

1 


i + 10.56 + 17.68 + 70.9 + I 


1 + 10.56 + 17.68 +70.9+0 2 

For this root, the second approximation was determined by dividing the 
coefficient E by the fourth sum from the quotient 


- 9.10 

70.9 

This procedure generally provides a good second approximation for the 
small real root. 

The cubic equation obtained by setting 

x 3 + 10.56X 2 + 17.68X + 70.9 

equal to zero is of the form 

-1- cX + d = 0 

In most lateral-stability work, a real root of this equation will be 
approximately equal to -b or it may be more closely approximated by 
the equation 


X = 


+ d 
+ c 


or for the particular case 


= _ ( 10. 56) 3 + 7 0. S . = 6 

(10.56) 2 + 17.68 95 



NACA TN 2409 


63 


Approximating the root by synthetic division 


1 + 10.56 + 17.68 + 70.9 

Approximation 

- 9.48 - 10.20 - 70.9 

-9 . 485 6 

- 9.49 - 10.16 - 71-4 

-9.49 5 

- 9.48 - 10.25 - 70.4 

-9.48 4 

- 9.45 - 10.50 - 67.9 

-9.45 3 

- 9-55 - 9-64 - 76.8 

-9.55 2 

- 9.65 - 8.78 - 85.9 

-9.65 1 

l + 0.91 + 8.90 - 15.0 

1 

l + l.Ol + 8.04 - 5-9 

2 

l + l.ll + 7.18 + 3.0 

3 

1 + 1.08 + 7.43 + 0.5 

4 

1 + 1.07 + 7.52 - 0.5 

5 

l + 1.075 + 7-48 0 

6 


For this large real root there is no simple method of determining the 
second approximation as there was in the case of the smaller real root. 
The magnitude of the estimated root in this case is arbitrarily 
increased or decreased slightly from the first approximation. From the 
remainders determined from the first two approximations, a fairly close 
third approximation can then be made . 

Factoring the quadratic equation obtained by setting 

X 2 + 1.0754 + 7.48 


equal to zero by use of the quadratic formula gives the final two roots 
of the biquadratic equation. ' ; 


X = - 


1.075 ±n/Fi 6 - 29.92 


= -0.538 


+ . fie: 

± ^-4 


76 


= -0.538 ± 2.68i 
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The roots of the biquadratic equation may be checked by multiplying 
the four factors to determine whether their product equals. the original 
biquadratic . , 

(X. - 0.1284)(x + 9.485)(X + 0,53.8 + 2.68i)(x + 0.538 - 2.68i) = (X 2 + 9.457X. 
1 . 220 ) ( X 2 + 1.07X + 7.47) = x 1 * + 10.43x3 + 16.32X 2 ■.+ 6.8. 6x - 9.10 

Solution by Descartes' Method 

Descartes' method of solving a biquadratic equation is particularly 
useful for solving equations which do not have any real roots, This 
method' is explained in most text books on advanced algebra and theory of 
equations. In general, the method consists of reducing the biquadratic 
equation to a'ciibic equation which can be solved easily. One root of 
the cubic equation is used to form two quadratic equations the roots of 
which are used to obtain the roots of the biquadratic equation. 


Method . - Reduce the general biquadratic equation 
AX^ + BX 3 + CX 2 + DX + E = 0 


tc the form 


+ bA,^ 4- c\^~ + d\ + e = 0 


by dividing by A. 

Obtain the values of q, r, and s from the following equations: 

q = c - 2 b 2 

8 

„ , be 1 . 3 

r.d -- + g b 



b 2 c 3 

" 25^ 



s 


e - 
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and form the equation 






O P / 

and solve this cubic equation in x for one of its roots f 0. 
Solution by trial by means of synthetic division is recommended. 
Determine the values of Z. and m from the equation 



1 

2 



r 


<1 0 2 r 

m = 2 + 2x + ^ 

Substitute the values of Z and m and the value of x used in 
obtaining Z and m in the equations 

2 , 

y + 2xy + Z = 0 ‘ 


■ ' . p 

y - 2xy -f m = 0 

and solve these quadratic equations for their roots \ y from which the 
roots of the biquadratic equation may be obtained from the following 
relation: 
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SPECIAL NOTATION USED IN CALCULATING MOTIONS WHEN 
THE CHARACTERISTIC EQUATION HAS COMPLEX ROOTS 


When two of the roots X]_ and X 2 are conjugate complex, the coefficients A]_ and Ag, 
Bi and B2> Ci and C 2 will he conjugate complex. If R + Ii is one of the roots Xi and 
if the powers of X]_ are expressed as 


then 



+ V 



Rl + l x i 

R 2 + ^2^ 
R3 + 131 

R 4 + * 4 * 


Substitution of the root R + Ii in the expression for Aq gives 


^qR^ + aqRj^ + a 2^3 a 3^2 a 4 R i + * 

> 

a 5 <1 

1 + ^IqIc^ + a^_l4 + ELryL^ + a3l2 + a||lqj 

1 i 

1 

'6AR 5 + 5BR4 + 4CR 3 + 3DR 2 + 2ER x J 

\ + i 

[6AI5 + 5BI4 + 4CI 3 + 3DI 2 + 2EI} j 

|i 


The division of these complex numbers is indicated by the equation 
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X 1 + ^l 1 _ x l x 2 + 3^2 *2^ ~ x 1^2 

x 2 + y 2 1 x 2 2 + yg 2 x 2 2 + *2 

It is evident from these relations that Ai is a complex number. In this case new symbols are 
used to represent the real and imaginary parts of A]_ as follows: 

A 1 = R A + I A i 

Ag is the conjugate of A 1 and will be referred to as 

A 2 = R A - I A i 

By procedures similar to those for the A coefficients. 


+ hjR]^ 4 * h^^ + + + 

b 5 ) 

+ 1 


ii 

1 

[6AR5 4 5 ERl. + + 31®2 + 2ER2] 

1 + 

(6a.I 5 + 5BI4 + 4 CI 3 + 3 DI 2 + 2 EI-J i 


which may he referred to as 


B 1 = R B + Igi 


and 


b 2 = r B " I B i 


Also, 


(cqR^ + c-^ + CgR^ + c 3 R 2 + c^R-l] 

1 + 

( C 0 I 5 + C 1 I 4 + C 2 I 3 + C 3 I 2 + c 4 I l' 

li 

(gAR^ + 5BR4 + 4CR 3 + 3DR 2 + 2ERJ 

1 + 1 

[6At 5 + 5BI4 + 4CI 3 + 3DI 2 + 2EI 1 j 

|i 


which may he referred to as 


o\ 
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and 


C'l - R c + I c i 


CT\ 

cx> 


c 2 “ R C " I C i 


Similar analysis shows that, if the roots and X .4 - are also conjugate complex quantities 

'^2 = R' + I'i and X .4 = R' - I'i), then 

a 3= r, A + i V 


and 


where 


\ “ R ’a - P A i 


^qR 1 ^ a-jR ' + a^R 1 ^ + a^R 1 ^ + . + a^) + + a^I'j^ agl* ^ + a^X^Ji 

(6aRV 5 + 5BR ! 4 +. 4CR*3 + 3BR|2 + 2ER’i) + (6ai ! 5 + 5BI r 4 + kci'^ + 3DI ! 2 + Sp^Ji ~ 


Also, 


and 


b 3 = R 'b + 1 'b 1 


B 4= R ’b - 1 'b 1 
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where 


^IdqR 1 ^ + "b^R * i|. + b 2 R ’ 3 + ^ 3^*2 b|^R 1 + 

b 5) 

+ 1 

[^q! 1 cj + ’bj_I , i| + b 2 I 1 3 + b^jl *2 + 

li 

i 

'6a R' 5 + 5BR’- 4 + 4CR'3 + 3DR*2 + 2ER']_j 

1 + 

(6AI'5 + 5BI'4 + 4CI'3 + 3DI'2 + 2EI'lJ 

|i 


Similarly, 


and 


where 


c 3 “ R 'c + I ’c i 
c 4 = R' C - I' c ± 


(c 0 R ! 5 + c l^'4 + c 2 ^ 1 3 + c 3^'2 + C 4 R r i) 

1 + 1 

^CqI 1 ^ + c-lI 1 !^ + c 2 I f 3 + c 3 I f 2 f c^I'J 

| i 

(6ar* 5 + 5 BR ' 4 + 4cr ' 3 + 3DR’ 2 + 2 ER , 1 ] + ^6ai’ 5 + 5 BI ' 4 + 4ci' 3 + 3DI' 2 + 2 EI'i] 



On 

VO 
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TABLE II.- REFERENCES CONTAINING USEFUL INFORMATION FOR ESTIMATING LATERAL STABILITY DERIVATIVES 
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Figure 2.- System of axes and angular relationship in flight. Arrows 
indicate positive direction of angles, t] = a - e. 
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Figure 4.- Lateral-stability boundaries calculated in table I. Ci 


the dependent variable. 


'*3 


was the independent variable. C 


3 


actually varied by changing Cy 


Varying C 
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in this manner 


was 


L Ptail' ~ ' U P 

caused changes in the tail contribution to all the other derivatives. 
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Figure 5*- Variation of lift-curve slope with aspect ratio, taper ratio, 
and sweepback for the case of subsonic incompressible flow. a Q = 0.11. 
Values from reference 3^- 
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Vertical location of horizontal tail, ZH/btai 

Figure 6.- Effect of horizontal -tail location on the effective aspect 

ratio of the vertical tail /A e , for the case of subsonic incom- 

V tail ) 

pressible flow, a = 0°. Taken from reference 35 * 





Aspect ratio, A 


Figure 7 -- Variation of with aspect ratio and sweep for the 

case of subsonic incompressible flow. X. = 1.0; § = 0. Taken from 

c 

reference 25. 



Aspect ratio , A 


Figure 8.- Variation of Cj / Cx with aspect ratio, taper ratio, and 

n 

sweep for the case of subsonic incompressible flow. Based on method 
of reference 25 . ' 
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Figure 9»- Effect of dihedral angle on for the case of subsonic 

incompressible flow. Taken from reference 5$ and 66. 
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Figure 10.- Variation of 


2 4 6 
Aspect ratio , A 
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ACL \ If CL \ with aspect 
, ^72 /\ °/a 


ratio for the case of subsonic incompressible flow, \ = 1,0. Taken 
from reference 96 . 
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Aspect ratio 


Figure 11.- Charts and formulas for estimating C 7 for the case of 

‘'D 


subsonic incompressible flow. Taken from reference 89 . 
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Theory (references 26, 34, and 107 ) 
Fairing based on 
experimental data 




Figure l4.~ Examples of suggested fairing of theoretical values of lift- 
curve slope for use in estimating values for the vertical tail in the 
transonic range. 
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